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Abstract 

We present a detailed study of charmless two-body B decays into final states involving two vec- 
tor mesons (1^1^) or two axial- vector mesons (AA) or one vector and one axial- vector meson (VA), 
within the framework of QCD factorization, where A is either a '^Pi or ^Pi axial- vector meson. 
The main results are as follows, (i) In the presence of NLO nonfactorizable corrections, effective 
Wilson coefficients are helicity dependent. For some penguin-dominated modes, the constructive 
(destructive) interference in the negative-helicity (longitudinal-helicity) amplitude of the B VV 
decay will render the former comparable to the latter and push up the transverse polarization, (ii) In 
QCD factorization, the transverse polarization fraction can be large for penguin-dominated charm- 
less VV modes by allowing for sizable penguin annihilation contributions, (hi) Using the measured 
K*^p~ channel as an input, we predict the branching ratios and polarization fractions for other 
B — > K*p decays, (iv) The smallness of the axial-vector decay constant of the ^Pi axial vector 
meson can be tested by measuring various bip modes to see if T{B^ — > p^) <^ T{B^ p ) 
and T{B^ P^) ^ T(B^ h\p^). (v) For the penguin-dominated modes aiK* and biK*, 

it is found that the former are dominated by transverse polarization amplitudes, whereas the lat- 
ter are governed by longitudinal polarization states, (vi) The rates of B ^ Ki{1270)K* and 
Ki{14:00)K* are generally very small. The decay modes KiK*~^ and Kj'K*~ are of particular 
interest as they are the only AV modes which receive contributions solely from weak annihila- 
tion, (vii) For tree-dominated B — > AA decays, the a^a^, a^a?, a^bf, a^b\, b^p~ and b\p~ 
modes have sizable branching ratios, of order (20 ~ 40) x 10~^. (viii) There are many penguin- 
dominated B AA decays within the reach of B factories: i^i(1270)(ai, 6^), ii'i(1400)(5i, a^), 
Ki(1270)(/i(1285),/i(1420)) and Ki(1400)(/i(1420), /ii(1170)). 
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I. INTRODUCTION 



Recently we have studied the charmless two-body B decays involving an axial-vector meson A 
and a pseudoscalar meson P in the final state [l|, Q]- There are two distinct types of axial- vector 
mesons, namely, ^Pi and ^Pi- We have studied their light-cone distribution amplitudes using the 
QCD sum rule method. Owing to the G-parity, the chiral-even two-parton light-cone distribution 
amplitudes of the ^Pi (^Pi) mesons are symmetric (antisymmetric) under the exchange of quark 
and anti-quark momentum fractions in the SU(3) limit. For chiral-odd light-cone distribution 
amplitudes, it is the other way around. In this work, we will generalize our previous study to 
charmless VA and AA modes. Moreover, we will use this chance to re-examine B — > VV decays. 

The charmless decays B — > VV, VA, AA are expected to have rich physics as they have three 
polarization states. Through polarization studies, these channels can shed light on the underlying 
helicity structure of the decay mechanism. Experimentally, B — > K*(p decays have been studied with 
full angular analysis and hence can provide information on polarization fractions and relative strong 
phases among various helicity amplitudes. Historically, it was the observation of large transverse 
polarization in i? — > K*(j) decays that had triggered a burst of theoretical and experimental interest 
in the study of charmless B — > VV decays. BaBar and Belle have observed that ~ 1/2 and /y ~ 
/_L ~ 1/4 in the K*(j) channels ^, Q], where fL,f± and /y are the longitudinal, perpendicular, and 
parallel polarization fractions, respectively. The transverse polarization fraction fx = + ~ 1/2 
is found to be of the same order magnitude as the longitudinal one /l in the penguin-dominated 
K*(p and K*p modes (except the decay B~ — > K*^ p^). While the naive expectation of /y ~ fj_ 
is borne out by experiment, the observed large Jt is in contradiction to the naive anticipation of 
a small transverse polarization of order /t ~ my/m\. This has promoted many to explore the 
possibility of new physics in penguin-dominated B VV decays. If so, the new physics effects 
should also manifest themselves in penguin-dominated VA and AA modes. 

The analysis of charmless B — > VV decays within the framework of QCD factorization 0] 
was first performed by us 0] followed by many others 0, fl Q Q Q, Q- In these studies, NLO 
corrections to the helicity-dependent coefficients such as vertex corrections, penguin contribu- 
tions and hard spectator scattering were calculated. However, most of the early results do not 
agree with each other due to the incorrect projection on the polarization states. Recently, Beneke, 
Rohrer and Yang 13| have used the correct light-cone projection operators and computed complete 
NLO corrections to and weak annihilation amplitudes. We will follow their work closely in the 
study of ^ VV decays. 

The generalization of the analysis oi B ^ VV decays to VA and AA modes is highly nontrivial. 
First of all, while the ^Pi meson behaves similarly to the vector meson, this is not the case for 
the ^Pi meson. For the latter, its decay constant vanishes in the SU(3) limit and its chiral-even 
two-parton light-cone distribution amplitude (LCDA) is anti-symmetric under the exchange of 
quark and anti-quark momentum fractions in the SU(3) limit due to the G parity, contrary to the 
symmetric behavior for the ^Pi meson. Second, there are two mixing effects for axial-vector mesons: 
one is the mixing between ^Pi and ^Pi states, e.g., Kia and Kib and the other is the mixing among 
^Pi or ^Pi states themselves. In this work we will derive the longitudinal and transverse projectors 
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for axial-vector mesons and work out the hard spectator scattering and annihilation contributions 
to VA and AA decays. 

Since the resolution of the K*(j) polarization anomaly may call for new physics beyond the 
standard model, this issue has received much attention in the past years. However, there are 
two crucial points that have been often overlooked in the literature. First, a reliable estimate of 
polarization fractions cannot be achieved unless the decay rate is correctly reproduced. Second, 
all the existing calculations except 0, @, [l^, 13] assume that the effective Wilson coefficients 
are helicity independent. This leads to the scaling law: fx ~ 0{my /m'^). Calculations based on 
naive factorization often predict too small B K*(j) and B K* p rates by a factor of 2 ~ 3. 
Obviously, it does not make sense at all to compare theory with experiment for fi^T at this stage 
as the definition of polarization fractions depends on the partial rate and hence the prediction 
can be easily off by a factor of 2 ~ 3. The first task is to have some mechanism to bring up 
the rates. While the QCD factorization and pQCD [l^ app roaches rely on penguin annihilation, 
soft-collinear effective theory invokes charming penguin [l5|] and the final-state interaction model 
considers final-state rescattering of intermediate charm states |l6l . llTl . Il8l | . Once the measured rate 
is reproduced, then it becomes sensible to ask what is the effect of this mechanism on polarization 
fractions. Next, it is important to consider NLO corrections to various helicity coefficients a^, 
such as vertex corrections, penguin and hard spectator scattering contributions. It turns out 
that in some of S — > VV decays, e.g. B K*cl), K*^ , NLO nonfactorizable corrections will 
render negative-helicity amplitude comparable to the longitudinal one and hence will bring up the 
transverse polarization. Therefore, any serious solution to the polarization puzzle should take into 
account NLO effects on of. 

There have been a few studies of charmless B — > AV and B AA decays in the literature 



19l . l2d . |21[ . Except for [19[] done in QCD factorization, the analysis in other two references was 



carried out in the framework of generalized factorization in which the nonfactorizable effects are 
described by the parameter A^^^, the effective number of colors. It has been claimed in [2l| that 
most of — > AV decays are suppressed and T{B — > AV) < T{B AP). This seems to be 
in contradiction to the naive anticipation that AV modes will have larger rates because of the 
existence of three polarization states for the vector meson. One of the main motivations for this 
work is to examine if the claim of 



211] holds. 



The present paper is organized as follows. In Sec. II we summarize all the input parameters 
relevant to the present work, such as the mixing angles, decay constants, form factors and light-cone 
distribution amplitudes for ^Pi and ^Pi axial-vector mesons and their Gegenbauer moments. We 
then apply QCD factorization in Sec. Ill to study B — > VV, VA, AA decays and derive the relevant 
spectator interaction and annihilation terms. Results and discussions are presented in Sec. IV. Sec. 
V contains our conclusions. Flavor operators and the factorizable amplitudes of selective B AV 
and AA decays are summarized in Appendices A and B, respectively. In Appendix C we give an 
explicit evaluation of the annihilation amplitude for the decay B VA. Since annihilation and 
hard spectator scattering amplitudes involve end-point divergences X^, we give explicit expressions 
of them for various VV, VA and AA modes in terms of X'^ in Appendices D and E. 
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II. INPUT PARAMETERS 



In this section we shall briefly discuss and summarize all the input parameters relevant to the 
present work, such as the mixing angles, decay constants, form factors and light-cone distribution 
amplitudes for vector and axial-vector mesons. 



A. Mixing angles 

Mixing angles of the axial- vector mesons have been discussed in [2^] and . Here we recapitulate 
the main points. For axial-vector mesons there are two mixing angles of interest: one is the mixing 
between ^Pi and ^Pi states, e.g., Kia and Kib and the other is the mixing among ^Pi or ^Pi states 
themselves, for example, the ^Pi states /i(1285) and /i(1420) have mixing due to SU(3) breaking 
effects. 

The non-strange axial vector mesons, for example, the neutral ai(1260) and &i(1235) cannot 
have mixing because of the opposite C-parities. In the isospin limit, charged ai(1260) and 61 (1235) 
also cannot have mixing because of the opposite G-parities. On the contrary, the strange partners 
of ai(1260) and 6i(1235), namely, Kia and Kib, respectively, are not mass eigenstates and they 
are mixed together due to the strange and non-strange light quark mass difference. We write 

i^i(1270) = KiAsinOKi + Kibcos9ki, 

ii:i(1400) = KiacosOki- KiBsmex,. (2.1) 



Various experimental information yields 9ki ~ ±37° and ±58° (see e.g. 23|). The sign of 9ki 
is intimately related to the relative phase of the Kia and Kib states. We choose the phase 
convention such that the decay constants of Kia and Kib are of the same sign, while the B — > Kia 
and B — > Kib form factors are opposite in sign. In this convention for Kia and Kib, the mixing 
angle is favored to be negative as implied by the experimental measurement of the ratio of Ki'j 
production in B decays [l|, 13] • 

Just like the rj — rj' mixing in the pseudoscalar sector, the l^Pi states /ii(1170) and /ii(1380) 
may be mixed in terms of the pure octet and singlet hi, 

|/ii(1170)) = \hi)coseip^ + |/i8)sin6lip^, |/ii(1380)) = sin6'ip^ + j/is) cos 6lip^ , (2.2) 

and likewise the l^Pi states /i(1285) and /i(1420) have mixing via 

|/i(1285)) = |/i)cos03p, ± |/8)sin03p^, |/i(1420)) = -|/i)sin03p, + |/8)cos03p, . (2.3) 



Using the Gell-Mann-Okubo mass formula J25l. |26|. we found that the mixing angles Oip^ and Osp^ 
depend on the angle and are given by []| 

Oip^ = 25.2°, Osp^ = 27.9°, for Ok^ = -37°, 

Oip^ ~ 0°, fcp^ = 53.2°, for Ok^ = -58°. (2.4) 
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B. Decay constants and form factors 



Decay constants of vector and axial-vector mesons are defined as 

(y(p,e)|g27/ift|0) = -i/vmye*, 
(^(^) Pi (p,e) 1 527/^7591 1 0) = i/3pj(ip^)m3p^(ipj)e*. (2.5) 
Transverse decay constants are defined via the tensor current by 

f«Pi(p,e)|^2a^.75gi|0) = fsi^p^^p" -elpn, 

{yip,^)\q2a,M0) = -f^ielp"" -elp^. (2.6) 
The decay constants /ip^ of the ^Pi non-strange neutral mesons 65(1235), /ii(1170), /ii(1380) vanish 
due to charge conjugation invariance. Likewise, the decay constant /^^ of the charged 6i vanishes 
owing to its even G-parity valid in the isospin limit. In general, the decay constants /ip^ and f^-p^ 
are zero in the SU(3) limit. As discussed in [it, they are related to f^p^ and /sp^, respectively, via 

where Oq are the zeroth Gegenbauer moment of <I>|| ^ to be defined later. Since we will assume 
isospin symmetry in practical calculations, this means that fip^ = for the bi and hi mesons and 
fi-p^ = for ai and fi mesons. Note that since /ip^ and fi-p^ are G-parity violating quantities, 
their signs have to be flipped from particle to antiparticle due to the G-parity, for example, fT^+ = 

—fj^- ■ In the present work, the G-parity violating parameters, e.g. af , af'^^'*, clq '^^^ ,cii'^^^ and 
^02^^ , are considered for mesons containing a strange quark. 

For the decay constants /jj(i285) •^/i(i420) I'^Pi states defined by 

(0k-7M75g|/i(1285)(P,A)) = -im^^(i285)/;,(i285)^^ (2-8) 
(0|w55l/i(1420)(P,A)) = -^m^l(l420)/;,(l420)^^ 
and the tensor decay constants for l^Pi states defined by 

{0\qcT^uq\hi{n70){P,X)) = ^4';'fn70) WefA)^^ (2-9) 
{0\qa^uq\hi{1380){P,X)) = ^/^^'(^aso) WfA)^^ 
the reader is referred to [l|, [^l for details. 

Form factors for the B ^ A and P — > F transitions read as 
{A{p,X)\A^\B{pb)) = i e,,^t,el^^^p%pfA^^{q^), 



niB - ruA 



_ r yba( 2\ 

{A{p,\)\V^\B{pb)) = -\{mB-mA)el,'>Vf^{q')-{e^^>-pB){pB+p)^^^ 

' ^ niB — rriA 



-2mA- 2^9^ 



{V{p,X)\V,\B{pb)) = -i e,,^pel^-^p%p^V"' {q% 

niB + Tny ^ ' 

{V{p,X)\A^\B{pb)) = (ms +mv)eW*yir(g2) _ (^W* . p^)(p^ + p)^/^''^'^') 



niB + mv 



2my-^q'^ Ar {q') - {q') , (2-10) 
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TABLE I: Form factors for B 
model 



ai,bi, Kiji, KiB transitions obtained in the covariant light-front 



27l | are fitted to the 3-parameter form Eq. (j2.12p except for the form factor V2 denoted by 



* for which the fit formula Eq. (j2.13p is used. 



F 


m 


a 


b 


F 


m 


a 


b 




0.25 


1.51 


0.64 




-0.10 


1.92 


1.62 


TrBai 
^0 


0.13 


1.71 


1.23 


yBbi 


-0.39 


1.41 


0.66 


j/Bai 
^1 


"^7 




14 


yBb^ 


18 


1 0"^ 


39 


T/Bai 
^2 


0.18 


1.14 


0.49 


yBb, 


0.03* 


2.13* 


2.39* 




0.26 


1.47 


0.59 


j\BK^B 


-0.11 


1.88 


1.53 




0.14 


1.62 


1.14 


^/BKiB 


-0.41 


1.40 


0.64 


yBK,A 


0.39 


0.21 


0.16 


yBK.B 


-0.19 


0.96 


0.30 


yBK,A 


0.17 


1.02 


0.45 


yBK,s 


0.05* 


1.78* 


2.12* 



where q = Pb — P, ^3^"^(0) 



yBA 



(0) and 



BAi 



2mA 

niB + my ^BV t 2\ 



Aril' 



ms + ruA ..BA 



2mA 
rriB — my 



(2.11) 



2my 2my 
Form factors for B ai(1260), &i(1235), Xi^, -fCiB transitions have been calculated in the rela- 
tivistic covariant li ght- front quark model (LFQM) (Table IH) [27|], the light-cone sum rule (LCSR) 
method (Table HI]) 28l |. and the pQCD approach [29|. Various B ^ A form factors also can be 



obatined in the Isgur-Scora-Grinstein-Wise (ISGW) model 30, l3l|] based on the nonrelativistic 
constituent quark picture. However, as pointed out in [l], the predicted form factor V^"'^{0) ~ 1.0 

^ afiT^ and the 



in the ISGW2 model 



3l[ | is too big and will lead to too large rates for B 
■0 



a^vr ) ^ B{B a-^ tt^), in contradiction to the experimental result 



wrong pattern B(B^ 

7r+). This may imply that relativistic effects in heavy-to-light 
transitions at maximum recoil that have been neglected in the ISGW model should be taken into 
account in order to get realistic form factors. 

It should be stressed that in the convention of the present work and LCSR, the decay constants of 
^Pi and ^Pi axial-vector mesons are of the same sign, while form factors V^'^ and have 
opposite signs. The sign convention is the other way around in the LFQM and pQCD calculations. 
Therefore, as explained in [ll, we put additional minus signs to the B ^Pi form factors in Table 

m 

The momentum dependence of the form factors calculated in the light-front quark model and 
the LCSR approach is parametrized in the three-parameter form: 

m 



F{q') 



(2.12) 



1 — aq^ /m\ + bq^/mj^ 

In the LFQM we use a different parametrization for the form factor ^2(5^) in some transitions [27] 

m 



F{q' 



(1 — 9^/"T'^)[l — aq^ /m\ + bq'^/m' 



(2.13) 



B\ 
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TABLE II: Same as Table U but in the light -cone sum rule model j28l |. 



F 


m 


a 


h 


F 


m 


a 


b 




0.30 ± 0.05 


1.64 


0.986 




-0.16 ±0.03 


1.69 


0.910 


j/Bai 


0.30 ± 0.05 


1.77 


0.926 


yBbl 


-0.39 ±0.07 


1.22 


0.426 


yBai 


0.60 ±0.11 


0.645 


0.250 


yBb^ 


-0.32 ±0.06 


0.748 


0.063 




0.26 ± 0.05 


1.48 


1.00 




-0.06 ±0.01 


0.539 


1.76 




0.27 ±0.05 


1.60 


0.974 


jyiPLiB 


^■^^-0.04 


1.72 


0.912 




0.22 ±0.04 


2.40 


1.78 


■trBKiB 


-o.45l°:J^ 


1.34 


0.690 


yBKiA 


0.56 ±0.11 


0.635 


0.211 


D jy 

yBKiB 


-o.48l°:S^ 


0.729 


0.074 


\/-BKiA 


0.25 ±0.05 


1.51 


1.18 


TrBKiB 
^2 


"-'•■'-^-0.02 


0.919 


0.855 




0.18 ±0.03 


1.63 


0.900 


j^Bhi 


-0.10 ±0.02 


1.54 


0.848 


yBji 


0.18 ±0.03 


1.81 


0.880 


j/Bhi 


-0.24 ±0.04 


1.16 


0.294 


R f 


0.37 ±0.07 


0.640 


0.153 


yBhl 


-0.21 ±0.04 


0.612 


0.078 




0.16 ±0.03 


1.47 


0.956 




-0.04 ±0.01 


0.500 


1.63 


ABfs 


U. lo It U.uZ 


1 (\A 


n Q1 Q 


ABhs 


— U.Uo It U.UZ 


1 

l.OD 


U.oZ ( 


yBh 


0.12 ±0.02 


1.84 


0.749 




-0.18 ±0.03 


1.22 


0.609 


yBf. 


0.26 ± 0.05 


0.644 


0.209 


yBhs 


-0.18 ±0.03 


0.623 


0.094 


yBh 
^2 


0.11 ±0.02 


1.49 


1.09 


TrBhs 

^2 


-0.03 ±0.01 


0.529 


1.53 



For B 
rules. 



p,K*,uj form factors, we shall use the results in 3^ obtained from light-cone sum 



C. Light-cone distribution amplitudes 

The light-cone distribution amplitudes (LCDAs) relevant for the present study are defined as 



22 



*(A) 



€ ^ 'Z 



1 e*Wz 1 



ly{P, A)|g-i(y)7^75g2(x)|0) = ifv f 1 - ^ ""^^ jl". ) W e(A>'^" 



fv mv 



1 (' 
4 



(2.14) 



(2.15) 



{V{P,X)\qi{y)a^,q2{x 



-fv / due 
Jo 



+ 



Hupy+up.) _ ££V/.)^±(^) 

mi e*(^)2; , -s , U) , ^ 

^^^y (PflZu - PvZ^) /l^i (U) 



p ■ z 



(2.16) 



V fv J Jo 



py+upx) 



hff\u) 



for the vector meson, and 
(^(P,A)|gi(y)7^ 75^2(2; 



{A{P,X)\qMlM^ 



(2.17) 



if Am A / dne^^^P^'+^P^) 

JO 

1 e*Wz 



2 '^(pz): 



-ifAmAe^upcT€*(X)P''z'' / dn 



/ rill p^iupy+upx) y± 



(2.18) 
(2.19) 



(.4(P,A)|g-i(y)a^,7592(x)|0) = /i / dtxe^("™+^P-) (el^^V.-^l'W)^^^ 



+ 



mi e*'^^h 



{pzf 



m. 



pz 



(^(P,A)|g-i(y)75g2(x)|0) = /imi(e*Wz) / dne*("™+^P^) 



(2.20) 



(2.21) 



for the axial-vector meson, where z = y — x with z"^ = and we have introduced the hght-hke 
vector p^ = — my^j^^^Zfj,/{2Pz) with the meson momentum P^ = rriy^^y Here $||, <I>j_ are twist- 

2 LCDAs, g^l\g^^\h\l\h^J'^ twist-3 ones, and 53, /13 twist-4. In the definitions of LCDAs, the 

longitudinal and transverse projections of polarization vectors e^^ along the z— direction for the 
(axial-)vector meson are given by [ssl] 



II M p2 



m 



ViA) 



2pz 



.*(A) _ *(A) _ *{A) 



(2.22) 



One should distinguish the above projectors from the exactly longitudinal and transverse polariza- 
tion vectors of the (axial-)vector meson, which are independent of the coordinate variable z = y — x, 
defined as 



.*(o)m 



E 



mv{A) _ 



m 



1- 



V{A) 



4^2 



^*(A)^ _± __ ^^^^^ 



(2.23) 

where we have defined two light-like vectors n^. with = (1,0,0,-1), and n'^ = (1,0,0,1) and 
assumed that the meson moves along the direction. 

In the QCDF calculation, the LCDAs of the vector meson appear in the following way [sS] 



{ViP,X)\qiaiy)q2sixm 



4 Jo 



du e 



i{u py+upx) 



8 



X { fymy i> 



pz 



51 W 



+0[{x-yf 



(2.24) 
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Here, all the components of the parton should be taken into account in the calculation before the 
collinear approximation is applied, so that one can assign the momenta 

^ ' ,n^t , = uEn^_ -k^^ + , (2.25) 



fcf = uEn^_ + A;^ + 



AuE ^ - X ' 

to the quark and antiquark, respectively, in an energetic light final-state meson with the momentum 
and mass m, satisfying the relation = En^ + m'^n^/{4:E) ~ En^. To obtain the light-cone 
projection operator of the meson in the momentum space, we apply the following substitution in 
the calculation 



d 



nt d 



dki. 



d 



2Edu dk^. 



(2.26) 



where terms of order k\ have been omitted. Moreover, to perform the calculation in the momentum 
space, we need to express Eq. ()2.24p in terms of z-independent variables, P and e*^'^\ instead of p 
and Consequently, the light-cone projection operator of the meson in the momentum space, 

including twist-3 two-parton distribution amplitudes, reads 

M5„ = M5„||+M5„^, 



(2.27) 



where Msa\\ and Msa± are the longitudinal and transverse projectors, respectively. 
For the vector meson, the longitudinal projector reads [34!] 



M, 



Jv "^v(e(A)"+) 



4 
iE 



$l,(u) 
Mi 



J^mv "^v(e(A)™+) 



dv{^^{v)-h\;'{v)) a^,n^^ 



(A)' 
2E 1 



dk 



+ 



+ 



k=up 



(2.28) 



and the transverse projector has the form 
Jvmv 



E 

'o 



dv{^\\{v)-g^l\v)) fc.ei; 



MX) 



d 



+8 4 dk^a 



dk^^ 
+ 



k=up 



(2.29) 



where k± is the transverse momentum of the qi quark in the vector meson. For the axial-vector 
meson, the longitudinal projector is given by 



-^T 2 ^ll(^) + 2E 



iE I dv {^±{v) - hf{v)) a^ulbW 



d 



dk< 



+ 75- 



+ o 



k=up 



,(2.30) 
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and the transverse projector given by 



Mi 



*(A) 



J Am A 



1 75 51 (w) 



E 



*(A) 







d 



dk 



I ■ u *(A)!^ p 



+ 



k=up 



(2.31) 



In the present study, we choose the coordinate systems in the Jackson convention; that is, in the 
B rest frame, one of the vector or axial-vector mesons is moving along the z axis of the coordinate 
system and the other along the —z axis, while the x axes of both daughter particles are parallel 



35| 



AO) 



{pc,0,0,Ei)/mi, e^^°) = (p^, 0, 0, -^2)^2, 
^M±i) ^ _L(o,^l,-i,0), 6^(±^) = -L(0,Tl,+i,0), 



*(±i) *(±i) 



where pc is the center mass momentum of the final state meson and • 
large energy limit, if the A meson moves along the direction, we will have e^^ 

*{A) 



*(A) 



(2.32) 

-5±i^±i. In the 
= 2EA/mA5x,o 



and • n„ =0. Note that if the coordinate systems are in the Jacob- Wick convention where 
the y axes of both decay particles are parallel, the transverse polarization vectors of the second 
meson will become e2(±l) = (0, ±1, — z, 0)/-v/2 and e*^^^^ • €2^^^^ = (5±i,±i. In general, the QCDF 
amplitudes can be reduced to the form of Jq dutr{M'^ . . .). 

To obtain the projector on the transverse polarization states in the helicity basis, one can insert 
^1 = ^° obtain 



Jvmv { 



(A)n ^ i i-")! \ _L 9T\'^) 

^1-75) gyw ± —. — 



+ 4^'^(i+75)ur(«)T 



51 [u) 



-E ^_(1 - 75) l^j^ dv{^\\{v) - g^l\v)) T 



i?^_(l+75) ( / dv{^\\{v) - g^l\v)) ± 



d 



dk±„ 
d 



+ 



k I =0 



(2.33) 



and 



4 + 



-75$l(w) 



9± ju) 



+ 4'^(1 + 75)UV" 



+E ^-(1-75) 



'dvi<^\\{v)-g^l\v))T 



5?(n)' 



d 



dk±^ 
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d 



dk_i_^ 



+ 



k I =0 



(2.34) 



Applying equations of motion to LCDAs, one can obtain the following Wandzura-Wilczek rela- 
tions in which twist-3 LCDAs are related to the twist-2 ones [331 via 



for vector mesons, and 

9± [u 



V 



U 



dv + 
- dv + u 



1 <S>ii(v) 



V 



dv 

V 

1 



+ ... , 



dv 



(2n- 1) 
2 



dv 



V 



V 

1 $ I 



V 



dv + u 



V 



^ 

« <I>ii('t; 



1 



+ ... , 

■dv + 
+ ... , 



u 



dv + 

) ^ 
- dv + u 



1 ^11 fv) 



+ ... 



w 



{2u - 1) 
2 



dv 



u 







dv 

V 

V 



+ ... , 



1 ^^{v) 



dv 



dv 



+ 



+ 



(2.35) 
(2.36) 

(2.37) 
(2.38) 

(2.39) 
(2.40) 
(2.41) 
(2.42) 



/o V 

for axial-vector mesons, where the ellipses denote additional contributions from three-particle dis- 
tribution amplitudes containing gluons and terms proportional to light quark masses, which we do 
not consider here and below. Eqs. (j2.35p - (|2.42p further give us 



4 



Jv 



1 $11 (tt) 



9± [v) (v) 



91 (V) 



u 

u 



du = ^>+(w) , 

du = -^^{v) , 



h't\v) = -2 



-du 



u 

it) I 



du{^±{u) - h\ '{u)) = vv 



u 

^ ^>_L(n) 
u 

u 



-du 



du 



du — V 



-2^v{v). 

1 $x(-») 
u 

1 $11 



(2.43) 



du 
-du 



u 



vv^y{v), 

--]^[v^-{v)-v^+{v)y. 



for vector mesons, and 



du = ^+{v) , 
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9± [V) (a) 



« $11 (u) 
u 



du = -<I>_(t;) , 



/i;fP)(,;) =_2 



-du 



u 



du{^±{u) - h\!\u)) = vv 



'dn($||(n)-5Sr)(^)) = i 



u 

u 



■du 
■du 
du — V 



1 ^±{u) 

U 

1 (-u) 



du 
du 



u 



vv^aiv), (2.44) 



for axial- vector mesons. 

After applying the Wandzura-Wilczek relations, the transverse helicity projectors (|2.33p and 
(j2.34p can be simplified to 



4 ' 
Jvmv J *(A), 



7^(1 ± 75) -n^ ^-(1±75; 



d 



+ 



fc I =0 



(2.45) 



and 



J Am A J *(A) 



M^(u) = iifE ^-75^1(n) 



4. ^+(^) 



75 



+4'^^-(^) 



7^(l±75)-nE ?^-(l±75: 



75 



+ 



I =0 



(2.46) 



From Eqs. (|2.28p - (|2.3ip and (|2.35p - (|2.44p . we see that $+ and $_ project onto transversely 
polarized vector or axial-vector mesons in which quark and antiquark flips helicity, respectively, 
while $„(a) projects onto longitudinally polarized vector (axial-vector) mesons in which either the 
quark or antiquark flips helicity. 

We next specify the LCDAs for vector and axial-vector mesons. The general expressions of 
LCDAs are 



$y(x,/x) = 6x(l — x) 



l + ^ar(^)C^/2(2x-l) 



n=l 



and 



^y{x,fl) = 3 



2x-l + £a^'^(/x)P„+i(2x-l) 



n=l 



(2.47) 



(2.48) 



for the vector meson, where Pn{x) are the Legendre polynomials. The normalization of LCDAs is 



/ dx^vix) = 1, / dx^jj{x) = 0. 
Jo Jo 



(2.49) 
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The explicit expressions of the LCDAs of axial-vector mesons have been discussed in details in 
We use 



5(2n -1)2-1 



'Pi. 



iPi 3 



5(2m -1)2-1 



for ^Pi mesons, and 



n=l 



3Pi 3 



(2.50) 



$11^1 (n) = 6uu\l + 3al' ^(2u-l) + a^ 



5(2u- 1)^ - 1 



(u) = 6uu \ ' + 30^'""^^ {2u - 1) + a^'"^^ 



3p, 3 



5(2u -1)2-1 



4. ^i(2n - 1) + 5: ay^^it,)Pn+ii2u - 1) 



n=l 



(2.51) 



for mesons. The normalization conditions are 

,'Pi 
''0 ' 



/ (ix<l>||^^ (x) = a, 

/ dx^f'{x) = l, 
Jo 

f dx^'/'{x)= [ dx<^>l^'{x) = 0. 
Jo Jo 



(ix^>|| i(x) = 1, 



dx<^l'''{x) = ai'''\ 



(2.52) 



It should be stressed that the LCDAs 'I'li^^ and <I>j^^ are defined with the decay constants 
f^p^ and /sp, , respectively, even though their corresponding normalizations are /ip^ and /s^^- As 
stressed in [l|, if we employ the decay constants fip and fih to define the the LCDAs <I>||'^^ and 



they will have the form 

^l^'in) = /ip,6txu|l + /.ip,^af''^^cf'(2n-l)|. 



fs\6nuf^l + M3P, 4''''cf\2u - 1)| 



(2.53) 



where /xip^ = l/ajj' ^ and //sp^ = I/oq' ^ which become infinite in the SU(3) limit. Therefore, 
it is most convenient to use Eq. ()2.50p for the LCDA <I)y^^ and (I2.51|) for the LCDA which 
amount to treating the decay constant of ^Pi as f^p^ and the tensor decay constant of ^Pi as /sp^. 
Of course, this does not mean that /ip^ {fi~p^) is equal to f^p^ {hp^)- 
For the B meson, we use the light-cone projector [34 1 



M, 



B ifsmB j 1 + ' 



I3a 



75 



(2.54) 



a/3 



The integral of the B meson wave function is parameterized as | 

'y^*f(p)^^, (2.55) 

I- p Ap 

where 1 — p is the momentum fraction carried by the light spectator quark in the B meson. 
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D. A summary of input quantities 



It is useful to summarize all the input quantities we have used in this work. 

For the CKM matrix elements, we use the Wolfenstein parameters A = 0.807 it 0.018, A = 
0.2265 ± 0.0008, p = 0.141+[{;[j^? and ?? = 0.343 ± 0.016 [36,]- The corresponding three unitarity 
angles are a = {90.7tii)°, = (21-7 ± 0.017)° andr = {67. 6tli)°- 



For the running quark masses we shall use [25|, [37 1 



mbinib) = 4.2 GeV, m6(2.1GeV) = 4.94 GeV, mfe(lGeV) = 6.34 GeV, 

mc{mb) = 0.91 GeV, ^^(2.1 GeV) = 1.06 GeV, mdl GeV) = 1.32 GeV, 

m,(2.1GeV) = 95MeV, m,(lGeV) = 118 MeV, 

md(2.1GeV) = 5.0 MeV, m„(2.1GeV) = 2.2 MeV. (2.56) 

Among the quarks, the strange quark gives the major theoretical uncertainty to the decay am- 
plitude. Hence, we will only consider the uncertainty in the strange quark mass given by 
ms(2.1GeV) = 95 it 20 MeV. Notice that for the one-loop penguin contribution, the relevant 



quark mass is the pole mass rather than the current one 38|]. Since the penguin loop correction is 
governed by the ratio of the pole masses squared Si = {fnf°^'^/mf"^'^)'^ [see Eqs. ()3.18p and (j3.20p 
below] and since the pole mass is meaningful only for heavy quarks, we only need to consider the 
ratio of c and b quark pole masses given by Sc = (0.3)^. 

The strong coupling constants employed in the present work are 

Q^(4.2 GeV) = 0.221, a^(2.1 GeV) = 0.293, 0^(1.45 GeV) = 0.359, a^(l GeV) = 0.495 . 

(2.57) 

For longitudinal and transverse decay constants of the vector mesons, we use (in units of MeV) 

/p = 216±3, /^ = 187±5, /;^. = 220±5, ^ = 215 ± 5, 

/^^ = 165±9, /i = 151±9, /i. = 185±10, = 186 ± 9 , (2.58) 



where the values of fv and (IGeV) are taken from 39t |. 

The decay constants /sp^ for oi , /i, /§ ^ and /f^^(lGeV) for bi, hi, hg obtained from QCD 



sum rule methods are listed in 



22l |. For the decay constants of Kia and Kib we use 



= 250 ± 13 MeV, Jk.^ = aj''^^^/^^ ^ -28 MeV, (2.59) 

where uses of f^ig — ^ 10 MeV [3] and the value of ajj'^^^ from Table IIIII have been made. 
Therefore, 

/i^i(i270) = -184 ± 11 MeV , /i^,(i4oo) = 177 ± 12 MeV , for Ok, = -37°, 
/i^i(i270) = -234 ± 15 MeV, fK,(im = 100 ± 12 MeV , for Ok, = -58°. (2.60) 



^ Recall that /g and /i are SU(3)-octet and -singlet states. 
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TABLE III: Gegenbauer moments of $_l and ^>|| for I'^Pi and l^Pi mesons, respectively, where 



"'^^^ and Oq^^^ are updated from the B — > Ki'y analysis and a\ 



,KiA „-L,KiA J,KiB 
) ""2 ' "2 ' 



and 

are then obtained from Eq. (141) in [24]. The scale dependence of Gegenbauer moments is 



referred to Eq. ([237]) . 





||,ai(1260) 
"'2 


an ^ 
"2 


3 p 

II f 1 

an ^ 
"2 


"2 




i Kje\ 


n no _i_ n no 
— U.U2 ih U.Uz 


— U.U4 zh U.Uo 


— U.U( ± U.U4 


—U.Uo ± U.Uo 


-0.30^°:^^ 


z.z oe V 


— w.ui in U.ul 


— U.UO in u.uz 


n -1- n iT^ 

— U.UO in U.UO 


n DzL -1- n 09 

— u.u^ in U.U-^ 






±,ai(1260) 


"1 


n '•'8 




±,KlA 

"0 


-L.ii'iA 
«2 


i LieV 


— 1.U4 dz U.o4 


— l.Uo =b U.oo 


— l.li dz U.oi 


— l.Uo zb U.4o 


26+°-°^ 
'-'■^"-0.22 


0.02 ±0.21 


9 9 PpV 
z.z oe V 


— W.Ol HZ W.ZiU 


— U.OZr ZIZ W.ZiO 


n Rfi -U n 9zL 

— W.OU ZIZ W.Z('4: 


— vj.O'4: ZIZ yj.oi 


^•^"^-0.21 


0.01 ±0.15 




11,61(1235) 
"1 




11,4''' 

"1 


"1 


"0 


aP- 


1 GeV 


— 1.95 ± 0.35 


— 2.00 ± 0.35 


— 1.95 ± 0.35 


— 1.95 ± 0.45 


-0.15 ±0.15 


o.o9l°:li 


2.2 GeV 


-1.56 ±0.28 


-1.60 ±0.28 


-1.56 ±0.28 


-1.56 ±0.36 


-0.15 ±0.15 


0.06^°:}^ 




±,61(1235) 


^2 


02 * 






1 GeV 


0.03 ±0.19 


0.18 ±0.22 


0.14 ±0.22 


-0.02 ± 0.22 


0.30^°:°? 


2.2 GeV 


0.02 ±0.14 


0.14 ±0.17 


0.11 ±0.17 


-0.02 ±0.17 


0.251H 



The Gegenbauer moments a-"*"^'^ and a| '^"'"•'''^ have been studied using th e Q CD sum rule method. 



Here we employ the most recent updated values evaluated at /i = 1 GeV [40] 
af * = 0.03 ± 0.02, a^'^* = 0.04 ± 0.03, af * = 0.11 ± 0.09, a^'^* 



0.15 ±0.07, a2 



0.14 ±0.06, a? = 0.18 ±0.08, 



= 0.10 ±0.08, 
0.14 ±0.07. (2.61) 



Note that aX = 0, a:j'"'^=0 for V = p,u}, (j). The Gegenbauer moments a| ^'^^'^ for axial-vector 
mesons are ;u,n,nari;ed in Table m This table is taken from Q with some updates on the 



Gegenbauer moments a^'^^"^, ajj'^^^, a|' 



KiA 



,KiA 



, u,2 



,«2 



^-L.-ft'iB ^ Qjj'-.ai ^j^^ a^' ^ .As stressed 



before, the values of the G-parity violating Gegenbauer moments (e.g. af , af'^''*, a. 



,KiA 



and aQ2^^) are displayed for the mesons containing a strange quark. Their signs are flipped for 
the mesons containing a s quark. In general, $_ft:(l — x) = $^(x). 
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For the B meson, we shall use Ab(1 GeV) = (250 it 100) MeV for its wave function and /b = 
210 ± 20 MeV for its decay constant. 

The Wilson coefficients Cj(/i) at various scales, // = 4.4 GeV, 2.1 GeV, 1.45 GeV and 1 GeV 
are taken from [4l|. For the renormalization scale of the decay amplitude, we choose jj, = mh{mb). 
However, as will be discussed below, the hard spectator and annihilation contributions will be 
evaluated at the hard-collinear scale Hh = Vl^^h with ~ 500 MeV. 



III. B VA, AA DECAYS IN QCD FACTORIZATION 



Within the framework of QCD factorization 
written in the form 



the effective Hamiltonian matrix elements are 



V ^ p=u,c 



h,P\B) 



(3.1) 



where Xp = VpbV*g with q = d,s, and the superscript h denotes the helicity of the final-state meson. 
Tj^'^ describes contributions from naive factorization, vertex corrections, penguin contractions and 
spectator scattering expressed in terms of the flavor operators af''^, while 7g contains annihilation 
topology amplitudes characterized by the annihilation operators 

The flavor operators a^'^ are basically the Wilson coefficients in conjunction with short-distance 
nonfactorizable corrections such as vertex corrections and hard spectator interactions. In general, 
they have the expressions @] 



p,h 



{M1M2 



Ci + 



+ 



Ci±l 

Ci±i Cpas 



iVf (M2) / ^^'^^^{x)dx 
Jo 

47^2 



(3.2) 



where i = 1, 



Nc 47r 

, 10, the upper (lower) signs apply when i is odd (even), Cj are the Wilson coefficients. 



Cf = {N^-l)/{2Nc) with = 3, M2 is the emitted meson and Mi shares the same spectator quark 
with the B meson. The quantities V/^{M2) account for vertex corrections, H-^{MiM2) for hard 
spectator interactions with a hard gluon exchange between the emitted meson and the spectator 
quark of the B meson and Pj(M2) for penguin contractions. The LCDA ^^^'^ in the first term of 
Eq. (|3.2p is for h = and for h = ±. The expression of the quantities NI^{M2) reads 

0, i = 6,8, 

1, else. 



iV/'(M2 



(3.3) 



Vertex corrections 



The vertex corrections are given by 



dx$;"2(x) 121n— -18 + 5(x) 
" L ^ 

dx^^^^ix) [-121n— + 6-5(1 -x) 
" L M 



y,°(M2) = 



dx^rn2ix) — 6 + h{x 



{i = 1-4, 9, 10) 
(i = 5,7) 
(i = 6,8) 



(3.4) 
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dx ^f-" (x) 12 In ^ - 18 + gT{x) 
L M 



= 1-4, 9, 10) 



dx^^^{x) -121n — + 6-5r(l-ic) , (« = 5,7) 

L ^ J 



0, 



(i = 6,! 



(3.5) 



with 



/ 1 - 2x , 
g\x) = 3 In X — ITT 



+ 



1 — X 

2Li2(x) -ln2a;+ 

1 — X 



(3 + 2ot) ln.x-{x-^l-x) 



h{x) = 2Li2(a;) — In^x — (1 + 2m) In a; — (x 1 — x) , 

Inx 

ffTla;) = 9{x) + — , 



(3.6) 



where x = 1 — x, is a twist-2 light-cone distribution ampUtudc of the meson M, (for 
the longitudinal component), and ^± (for transverse components) are twist-3 ones. Specifically, 
$^ = for M = y and <^rn = for M = A. 

Hard spectator terms 

H^{MiM2) arise from hard spectator interactions with a hard gluon exchange between the 
emitted meson and the spectator quark of the B meson. Hf{MiM2) have the expressions: 



yiBMi,M2) As Jo 
^0 



uv 



uv 



(3.7) 



for i = 1 - 4, 9, 10, 



^(SMi,M2) As Jo 



(3.8) 



for i = 5,7, and Hf{MiM2) = for i = 6,8, where the upper (lower) signs apply when Mi = V 
(Ml = ^). The transverse hard spectator terms Hf{MiM2) read 



(3.9) 



ii;^(MiM2) = .^(MiM^)^^^^^^^^;^-"^^^"^^^- ^ /^.d.i^^^^l^£M<!M,(3.10) 



for i = 1 - 4, 9, 10, and 



(3.11) 



//+(MiM2) = a+(MiM2)'^^^^^-^^-"^^-^^-!?^ dud. i^i^^^li^^^^f^, (3.12) 



2 ^{BM^Mi) As io 
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for i = 5,7, and 

Hr{MiM2) 

H+{MiM2} 
for i = 6, 8, where 



a'{MiM2 



0, 



ifsfMifMi'mMi mBruMi ms 



ViIbX, 



M2 



A 



B Jo 



dudv 



vuv 



(3.13) 
(3.14) 



r]-{MiM2) 
a~{MiM2) 
a+{MiM2) 



+1 


for M1M2 


= VV, VA 


-1 


for M1M2 


= AV, A A, 


+1 


for M1M2 


= VA, AV 


-1 


for M1M2 


= VV, AA, 


+1 


for M1M2 


= VA, AA 


-1 


for M1M2 


= VV, AV 



(3.15) 



and r]~^ {M1M2) = — 1. To write down Eq. (|3.13p . we have factored out the r^'^^ term so that will 
contribute to the decay amplitude in the product of r^^^af' oc r^'^H^ . Two remarks are in order: 
(i) We have checked explicitly that the hard spectator terms depend on the B meson wave function 
$f (/?), but not on $f (p). (ii) Since Beneke et al. [13| adopted the Jacob convention for transverse 



polarization states, they have e* -62 = 1- As a consequence, their expressions for the parameters 77^, 
o"^ and the decay amplitudes x'^'^^'^'^^ defined below have signs opposite to ours. Nevertheless, 
the expressions for H^{MiM2) are independent of the choice for transverse polarization vectors. 
The helicity dependent factorizable amplitudes defined by 



xiBM,,M,) ^ (^M2{p2,e*2)\J^\0){AMp,,el)\J>^\B) 



(3.16) 



have the expressions 



{BVi,M2) 



X, 



(BAi,M2) 




2mvi 



(m| 



(m| 



'^M2)i^B + mvi)^f^^(g^ 



BVi („2\ ^ '^Pc 



X. 



(BVi,M2) 



X 



(BAi,M2) 



[(- 




rriB ) 


[(- 






ruB ) 



rtiB + 
2pc 



rriB + 

rriB - 
V'^'^Hq^) 



rriB - ruAi 



(3.17) 



-ifMiiTT-BrriM. 

-ifMiirLBrnM. 
where M2 = V2, A2 and Pc is the cm. momentum. 

Penguin terms 

At order a^, corrections from penguin contractions are present only for i = 4, 6. For i = 4 we 
obtain 



P4'^(M2) 



Cl 



{sp) + + C3 Gl. (ss) + G\, (1) + 2g 



+{cA + ce)J2 [Gii^i^i) + 9'm2 



(3.18) 
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where Si = mf/ml and the function G\r{s) is given by 



9M2 





4 rub 2 

- in 1 — 

3 /X 3 



dx XX ln[s — uxx — ie] , 
^^^'''{x)dx, 



dm 



4 mfe 
- in — 
3 fi Jo 



<^^^^'^{x)dx, 

with $^^2,0 ^ (^Ma^ $^2,± = ^M2_ i = 6, the result for the penguin contribution is 



^6 



(M2) 



A-hiSp, 



+ C3 



^M2 



+ {ci + Cf,)Y,Gli^{si 



(3.19) 



(3.20) 



In analogy with (|3.19|) . the function Gm^{s) is defined as 



4 
0. 



du ^>m2 (n) 



ln[s 



uxx 



(3.21) 

Therefore, the transverse penguin contractions vanish for z = 6, 8: Pq^^ = 0. Note that we have 
factored out the term in Eq. (j3.20p so that when the vertex correction Vg^g is neglected, Og 
will contribute to the decay amplitude in the product r^'^a^ 



pO 



For i = 8, 10 we find 



-'^10 




{ci+NcC2)Gii^{Sp,, 

Gii^{sp) + 2gM2 



For i = 7, 9, 



(M2 



Sir 



m 



M2 



27tt 



(ci +iVcC2) 



5pc In ^ + (5p„ In ^ + 1 



(3.22) 
(3.23) 

(3.24) 



if M2 = p^,ll!,(P, otherwise P^^{M2) = 0. Here the first term is an electromagnetic penguin 
contribution to the transverse helicity amplitude enhanced by a factor of mBmh/mj^^, as first 
pointed out in 4^. Note that the quark loop contains an ultraviolet divergence for both transverse 
and longitudinal components which must be subtracted in accordance with the scheme used to 
define the Wilson coefficients. The scale and scheme dependence after subtraction is required to 
cancel the scale and scheme dependence of the electroweak penguin coefficients. Therefore, the 
scale ^ in the above equation is the same as the one appearing in the expressions for the penguin 
corrections, e.g. Eq. (j3.19p . On the other hand, the scale z/ is referred to the scale of the decay 
constant /a/j (i^) as the operator q^y^q has a non- vanishing anomalous dimension in the presence of 
electromagnetic interactions 0]. The v dependence of Eq. (|3.24|) is compensated by that of fM2{^)- 
The relevant integrals for the dipole operators Og^-y are 



G^ 



1 ^t^'iu) 
du 





1 du 
u 



u 



M2I 



u^^^iu) + 



u^^\u)+u^'^''{u 



M2I 



(3.25) 
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Using Eq. ^M), Gf can be further reduced to 



G- 



du 



M21 



0. (3.26) 

Hence, G^ in Eq. (j3.25p are actually equal to zero. It was first pointed out by Kagan [ll!] that the 
dipole operators Q^g and Qj^ do not contribute to the transverse penguin amplitudes at 0{as) due 
to angular momentum conservation. 

Annihilation topologies 



The weak annihilation contributions to the decay B M1M2 can be described in terms of the 

(3.27) 



building blocks and ^^^^v 



V ^ p=u,C V ^ p=U,C i 



The building blocks have the expressions 



h 

^3,EW 
^4,EW 



h Gf r 



Gp j C'iT' 



Cf 



C2A\, 64 
Cg^i + C7(4 + ^3) + NcCfiA^ 



(3.28) 



where for simplicity we have omitted the superscripts p and h in above expressions. The subscripts 
1,2,3 of Aii^ denote the annihilation amplitudes induced from {V - A){V - A), {V - A){V + A) 
and (5 — P){S + P) operators, respectively, and the superscripts i and / refer to gluon emission 
from the initial and final-state quarks, respectively. Following 0] we choose the convention that 
Ml contains an antiquark from the weak vertex and M2 contains a quark from the weak vertex. 
The explicit expressions of weak annihilation amplitudes are: 

"1 



4'°(MiM2) 



1 



+ 



1 



u[\ — uv) ' uv^ 

2 ' 

uv 



A{-{MiM2) 
4' + (MiM2) 



-IT a. 



2mM^mM2 



dudvl ^^_!'{u)¥^-'{v) 



M21 



-TTQ 



2mM^mM2 



nil 



j\udv\^X'{u) ^X\v) 



U + V 



+ 



uv'^ {\ — uv)"^ v'^{l — uv) 



(3.29) 
(3.30) 
,(3.31) 



4°(MiM2) 



vra. 1^ dudv I $J^i(^;)$J^2(^) 



+ 



v{l — uv) u^v 



20 



2 

uv 



JO 

U + V 1 

+ 



^2 (M1M2) = -vra 



m% JO 
2 u 



(3.32) 



(3.33) 



(3.34) 



4'°(MiM2) 



2u 



+D 



uv{l — uv) 
2v 



4'-(MiM2) 



Jo 



^ II nv(l — uv) J 

I "T-Mi — UV) 



4'°(MiM2) = Tias j^dudvlc^''^^^^ r^^^ ^mM^^\''i.'^ 



uv{\ — uv) 
2{l + v) 



uv^ 



4' (M1M2) = -vra 



Jo \ ^ ^ 'u^v 



+ 



ruMi 



j^NHM, ^^M, ^M,(^) ^M.(^) 2 



and Af''^ = ^^''^ = 4+ 



: ^3'^ = 0, where 



^1 ' 


^3 • 




= D^^ = 


-1, 


^2 1 


^2 • 




= C^^ = 


-1, 


^3 


' ^3 1 ^3 


^AV 


= C^^ = 


-1, D^^ = D^^ 



(3.35) 



(3.36) 



(3.37) 



(3.38) 



-1, (3.39) 

and the parameters C and D are equal to +1 for all other cases. Note that our results for J^n^'^ 
have opposite signs to that in [13] as Beneke et al. adopted the Jacob convention for the transverse 
polarization vectors. We employ the same convention as in [6] that Mi contains an antiquark from 
the weak vertex with longitudinal fraction y, while M2 contains a quark from the weak vertex with 
momentum fraction x. 
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Since the annihilation contributions suppressed by a factor of mim2/m\ relative to 

f 

other terms, in numerical analysis we will consider only the annihilation contributions due to , 

^3 ' ^1,2,3 ^1^*-^ ^3 • 

Finally, two remarks are in order: (i) Although the parameters ai{i ^ 6,8) and a-Q^^ry- are 
formally renormalization scale and 75 scheme independent, in practice there exists some residual 
scale dependence in aj(/i) to finite order. To be specific, we shall evaluate the vertex corrections 
to the decay amplitude at the scale jj, = rrn,. (The issue with the renormalization scale /U will 
be discussed in more detail in Sec. IV). In contrast, as stressed in [S], the hard spectator and 
annihilation contributions should be evaluated at the hard-collinear scale = V /xA/^ with Kh ~ 
500 MeV. (ii) Power corrections in QCDF always involve troublesome endpoint divergences. For 
example, the annihilation amplitude has endpoint divergences even at twist-2 level and the hard 
spectator scattering diagram at twist-3 order is power suppressed and posses soft and collinear 
divergences arising from the soft spectator quark. Since the treatment of endpoint divergences is 
model dependent, subleading power corrections generally can be studied only in a phenomenological 
way. We shall follow 0] to model the endpoint divergence X = J^dx/x in the annihilation and 
hard spectator scattering diagrams as 

Xa = In (^^) (1 + pAc'^^), Xh = In (^^) (1 + pHe^'^^, (3-40) 

with the unknown real parameters PA,H and (j)A,H- For simplicity, we shall assume that X\ and 
X^ are helicity independent; that is, X^ = X^ = X\ and Xjj = X'^ = X^j. 



IV. NUMERICAL RESULTS 

The decay amplitude of S — > M1M2 with M = V,A has the general expression of 
^Mi('^Afi)eM2(AM2)-^/xi/ with Xmi,M2 being the corresponding helicities. Hence, the decay ampli- 
tude can be decomposed into three components, one for each helicity of the final state: Ao,A+,A-. 
The transverse amplitudes defined in the transversity basis are related to the helicity ones via 

A+ + A- A+-A- 



The decay rate can be expressed in terms of these amplitudes as 

r = /VdAP + \A+\' + \A.\') = + |-4||r + \A^n, (4.2) 

with pc being the cm. momentum of the final-state meson. Polarization fractions are defined as 

r \A P 

r |AP + |^|||2 + |^x|2' ^^-""^ 

with a = L, II, _L. The relative phases are 

^± = aTg{A±/Ao), (pi\ = arg(^||/A)- (4.4) 

Note that the experimental results of and (p± obtained by BaBar and Belle are for B (pK* 



decays [49|, l50|, ISll- According to the convention given by BaBar and Belle, |^+| > |^_| and 
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't'W = (j)± = IT for B — > <j)K* in the absence of final-state interactions. Since our calculations are 
for B (j)K* decays, in Eq. (I4.14|) below we shall transform BaBar and Belle results from 0|| to 
vr — (/)|| and (l)± to — (/>_l so that [.A+l < |^_| in i? ^ <pK* . When strong phases vanish, (/)|| = 0, 
(f)± = — vr for B (j)K* . 



A. B ^ VV decays 

The branching ratios and polarization fractions of charmless B VV decays have been mea- 
sured for pp, pco, pK*, (j)K*,ioK* and K*K* final states. It is naively expected that the helicity 
amplitudes Ah (helicities /i = 0, — , + ) for i? — > VV respect the hierarchy pattern 

Ao:A^:A^ = l:(^) : (^)\ (4.5) 
\ nib J \ rrib J 

Hence, they are dominated by the longitudinal polarization states and satisfy the scaling law. 



namely 11 1 



l-f,=0 ^], i±^l^o(^). (4.6) 

J f\\ KmBJ 

with fL,f± and /y being the longitudinal, perpendicular, and parallel polarization fractions, re- 
spectively. In sharp contrast to the pp case, the large fraction of transverse polarization observed 
in B ^ K*p and B K*(j) decays at B factories (see Table IIVI below) is thus a surprise and 
poses an interesting challenge for any theoretical interpretation. Therefore, in order to obtain a 
large transverse polarization m. B ^ K*p, K*(j), this scaling law must be circumvented in one way 



or another. Various mechanisms such as sizable penguin-induced annihilation contributions 111 ], 
final-state interactions 16, 18], form- factor tuning 43] and new physics 0, [l2, 44, 4^ (where only 
the models with large scalar or tensor coupling can explain the observation for f± ~ /y l^, |4j]) 
have been proposed for solving the B — > VV polarization puzzle. It has been shown that when 
the data for (l)K* and Krj^'^ modes are simultaneously taken in into account, the standard model 
predictions with weak annihilation corrections can explain the observation, while the new physics 



eff'ect due to scalar- type operators is negligible [46]. 



Before proceeding, we would like to make a few remarks on the polarization anomaly. First, 
the hierarchy of helicity amplitudes given by Eq. (j4.5p is valid only for factorizable VK-emission 
amplitudes. It may be violated in the presence of nonfactorizable corrections (e.g. vertex, penguin 
and hard spectator scattering contributions) and annihilation contributions. Indeed, we shall show 
below that the polarization pattern (|4.6|) will get modified when nonfactorizable contributions are 
included in QCD factorization. We shall see later that the polarization anomaly is not so serious 
as originally believed. Second, it is known that the predicted rates for the penguin dominated 
B — > VP, VV decays in QCD factorization are generally too small by a factor of 2 ~ 3 compared 
to the data. It is obvious that in order to have a reliable calculation for polarization fractions, 
it is of great importance to first reproduce the decay rates correctly. Otherwise, the estimation 
of fL,\\,± will not be trustworthy. Hence, our first priority is to have a mechanism resolving the 
branching ratio puzzle for the penguin dominated charmless B — > VV decays and hopefully the 
same mechanism also unravels the polarization anomaly. 
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1. Tree- dominated decays 



Branching ratios and polarization fractions for tree-dominated B ^ pp and puj are shown 
in Table IIV[ The theoretical errors correspond to the uncertainties due to variation of (i) the 
Gegenbauer moments, the decay constants, (ii) the heavy-to-light form factors and the strange quark 
mass, and (iii) the wave function of the B meson characterized by the parameter A^, the power 
corrections due to weak annihilation and hard spectator interactions described by the parameters 
PA,H, <t'A,H^ respectively. To obtain the errors shown in Tables IVTIllXIl we first scan randomly the 
points in the allowed ranges of the above nine parameters and then add errors in quadrature. More 
specifically, the second error in the table is referred to the uncertainties caused by the variation 
of PA,H and (pA.H, where all other uncertainties are lumped into the first error. Here we consider 
the default results for tree-dominated decays by setting the annihilation parameters to be zero, i.e. 
PA = (pA = 0, though the predictions are insensitive to the choice of them. 

It is obvious from Table HVl that the longitudinal amplitude dominates the tree-dominated decays 
except for the p^uj mode where the transverse polarization could be equally important. The naive 
expectation of ~ 1 — 4m^/m^ ~ 0.92 is experimentally confirmed. The calculated rates are also 
in agreement with experiment except that the predicted rate for B~ — > p~uj is slightly high. Its 
decay amplitude reads 



X 



(Bp,Lj) 



+ 



oP,h 



X 



{Bu,p) 
h 



(4.7) 



It is obvious that this decay is dominantly governed hy B ^ uo transition form factors. The data 
of B~ — > p~oj and B — > K*lv to be discussed below seem to suggest that B ^ uj form factors are 
slightly smaller than what are expected from the light-cone sum rules ^3|. 



2. Penguin-dominated decays 



The decays of interest in this category are B K*p, K*(j), K*u) and K*K*. 
B K*p 

We first consider B K*p decays. Retaining the leading contributions, their decay amplitudes 
are approximated by 

c,h 



^K'-pO 



V2Ab 

^^^K'-p+ 



Vua\ + V,{af + Hi) 



^pK'^ 



+ ^C2"3,EW 



X 



h 

K*p-> 



^K*p' 



(4i 



where Vp = VpbV*g with \ Vc\ ^ \Vu\, characterizes the penguin-induced weak annihilation (see 
Eq. dHU) for definition) and X^^, is a shorthand notation for Xg_^^j^, with its explicit expression 
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TABLE IV: CP-averaged branching ratios (in units of 10 ^) and polarization fractions for B 
pp, puj, K* p, K*(j), K*uj, K* K* decays. The annihilation parameters are specified to be pA = 0.78 
and 4>A = -43° for K*p,K*K* and pA = 0.65 and (t>A = -53° for K*(f) and K*uj by default. For 

and 



longitudinal polarization fraction, the theoretical uncertainty is dominated by pa,h and 4>a 
hence only this error is listed in the table fo r /r . Experimental results are taken from [4 



nence only 

3, 0, m i 



53|, 



54 



55 



JL 



47 



48 



5g, |57|, |58|, |59|, [GO, l6l| and the world averages from [62] . 



1 Icr* QTT 




B 






/l 








Theory 


Expt 


Theory 


Expt 


Theory 


Expt 


B- 




nn n+4.0+2.0 
^'-'•"-'-1.9-0.9 


18.2 ±3.0 


0.96tn^ 


n 01 9+0.044 


0.02 ±0.01 






p+p- 


oc; c+1.5+2.4 
^^•"3-2.6-1.5 


94 9+3-1 
^^•^-3.2 




n 078+0.025 
U.y (0_Q Q22 


"-'•"-'^-0.00 






p^p^ 


n q + 1.5+1.1 

U.y_o.4-0.2 


0.68 ±0.27 




0.70 ±0.15 






B- 


p~uj 


-IQ 9+3.3+1.7 
J^y-^-1.6-1.0 




0.96tn^ 


0.82 ±0.11 


0.02 ±0.01 






p^uj 


n 1+0.1+0.4 

u-^-o.i-o.o 


< 1.5 


U.OO_Q 29 




22+°-^^ 

"-'•^^-0.23 




B- 


K*^p- " 


Q 9+1.2+3.6 
^■^-1. 1-5.4 


9.2 ± 1.5 




0.48 ± 0.08 


'-'•^"3_o.26 




B- 


K*-p^ 


c; c+0.6+1.3 
^•"3-0.5-2.5 


< 6.1 




n qf;+0.06 b 
U-^"-0.16 


'-'•■'-O-0.15 






K*-p+ 


Q+1. 1+4.8 
0-^-1.0-5.5 


< 12 


^•'J'J-0.32 




'-'•^^-0.22 








A R+0. 6+3.5 
*-"-0. 5-3.5 


5.6 ± 1.6 


o.39i°:^; 


0.57 ±0.12 


U-'J'J-0.30 




B- 


K*-(l) ^ 


10 n+1-4+12.3 
^'-'•U-1.3- 6.1 


10.0 ± 1.1 


0.49^0;^^ 


0.50 ±0.05 


U.ZO_Q 25 


0.20 ±0.05 






+1.3+11.9 
^•"J-1.2- 5.9 


9.5 ±0.8 


IJ-OU.Q 42 


0.484 ± 0.034 


25+°-2l 

U.ZO_o.25 


0.256 ±0.032 


B- 


K*-uj 


q c;+0.4+3.0 
•-••^-0.4-1.7 


< 3.4 


'-'•"30_o.38 




1 7+°-20 








q n+0.5+2.9 
•J-U_o.4-1.8 


< 2.7 


"•'J' -0.46 




^•^J--0.22 




B 


K*^K*- 


n (^+0.1+0.3 
U-O-0. 1-0.3 


< 71 


'-'•*^-0.38 




n 07+0.19 






K*-K*+ 


n 1+0.0+0.1 

'-'•-'--0.0-0.1 


< 141 


1 













n f^+0.1+0.2 

"-'•"-0.1-0.3 


1 90+0.37 

-'-•^"-0.32 


^•"^^-0.48 


0.801°:}! 


"-'•^^-0.24 





"This mode is employed as an input for extracting the pararnaters pA and 4'a for B — > K* p decays. 
''A recent BaBar measurement gives fL{K*~p^) = 0.9 ± 0.2 [47|, but it has only 2.5(t significance. 
'^This mode is employed as an input for extracting the parameters pA and (j)A for B K*(j) decays. 



shown in Eq. (j3.17p . The expressions of the flavor parameters a^'^ in terms of the coefficients 
a^*'^ can be found in Eq. (lAlh . To proceed, we shall first neglect annihilation completely by 

h,p 



setting = 0. In the absence of NLO nonfactorizable corrections, the parameters 
independent and hence the hierarchy relation 



\X 



\x 



1 : 0.26 : 0.03 and IX^* 

I pK* 



are helicity 

3|) for helicity amplitudes is respected as 1-^^^*^! : 



\X 



\X 



pK* 



1 : 0.30 : 0.005. When vertex. 



p,h 



penguin and hard spectator corrections are taken into account, we see from Table IVl that a; 
and 03 EM ^^'^ negative helicity differ significantly from that the longitudinal ones. For example, 
the real parts of 02 and a^^^y^ have opposite signs for /i = and h = — . Let us consider two 
extreme cases for the longitudinal polarization fraction. From Eq. (14.81 ) we have 

3 , 



BO- 
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TABLE V: Longitudinal- and negative-helicity amplitude parameters. 
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(4.9) 



From Table IVl we see that the interference between 04" and Osevv constructive for h = — and 
destructive for /i = for the decay 



K*^p^ and the other way around for B 



K* 



P 



As a consequence, A" is comparable to A^ for the former but is highly suppressed relative to A^ 
for the latter. The longitudinal polarization fraction for the penguin dominated processes can be 
approximated as 

2 
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Q4' + CvO:^-^^ 



h=0, 
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h=0-|a3+"4 +2"3,EW + P: 
\2 aB^P ,t:^B 
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(4.10) 



rriK'/mBfAQ^P/F^-'P oc (mi^./mB)^ c„ = for and 



modes, c„ = 1 for i^T* p*^ and Cy = —1 for K*^ p^ (see Ref. 13|] for the definitions of the ^0 and 
form factors). The calculated branching ratios and the longitudinal polarization fractions fi in 
QCDF are shown in the case (i) of Table ED Indeed, we find fL{K*^p^) = 0.46 and fL{K*-p^) = 
0.97. If the coefficients are helicity independent, we will have fL{K*^ p^) = 0.91 rather than 
0.46 ! However, the NLO corrections to will render the negative helicity amplitude A~{K*^ p^) 
comparable to the longitudinal one AP{K*^p^) so that even at the short-distance level, for 
B^ K*^p^ can be as low as 50%. Similar detailed discussions for (j)K* modes will be given latter. 

Comparing with the data, it appears that even though the naive estimate of fi is too large for 
K*^p~, the experimental observation of a large fi for K*~ p^ and a small fL for K*^p^ are well 
accommodated. However, as stressed before, in order to have a trustworthy estimate of polarization 
fractions one has to reproduce the rates correctly as the predicted branching fractions for K*^p~ and 
K*^p^ are too small compared to experiment (see Table rvT|) . In the present work, we shall follow 



111 ] to ascribe the necessary enhancement to a potentially large penguin annihilation characterized 
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TABLE VI: CP-averaged branching ratios (in units of 10^^) and the longitudinal polarization 
fraction fi for B — > K*p and K*(f) decays for three cases: (i) no annihilation contribution, (ii) 
adding annihilation contributions with pA = 0.78, (j)A = —43° for K*p and pA = 0.65, (pA = —53° 
for K*(j). The predictions by Beneke, Rohrer and Yang [3] are shown in the last two columns. 
For simplicity, only the central values are exhibited. The theoretical uncertainties in case (ii) are 
shown in Table IIVI Experimental results are taken from 47 , 48'] . 
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by the parameter /33. We fit the data of B^ K*^ p^ by adjusting the parameters pA and (/>^ that 
characterize the nonperturbative effects of soft gluon exchanges in annihilation diagrams. Prom Fig. 
[T]^a) we see that pA is preferred to be around 0.78 in order to fit the rate, while the corresponding 
fi is around 0.48 (see Fig. [I^c)). Only the theoretical uncertainty due to the variation of the phase 
(j)A is considered in Fig. [TJ It is clear that the total branching ratio and the longitudinal one Bl 
increase with pA, whereas fi decreases slowly with pA- To fit the rate and /l simultaneously for 
B~ K*^p~, we find pA ~ 0.78 and (j)A ~ —43°. Using this set of parameters, we are able to 
predict branching ratios and polarization fractions for other B K* p decays as exhibited in Table 
IIVI and in case (ii) of Table IVll In the presence of penguin annihilation, the parameter ctj'* in Eqs. 
gSD and (gJD should be replaced by a/" + /3f . From Table |Vl one can check that both fL{K*^p^) 
and fL{K*^p^) will be decreased when penguin annihilation is turned on. 

Within the QCDF framework, Beneke, Rohrer and Yang (BRY) have employed the choice 
PA = 0.6 and 4>a = —40° obtained from a fit to the data oiK*<j) to study K* p decays [135 ■ They have 
noticed that the calculated K* p branching fractions are systematically below the measurements. 
This is not a surprise as their pA is smaller than 0.78 [see also Fig. (U^a)], since as emphasized before, 
the estimation of polarization fractions will not be reliable unless the calculated partial rate agrees 
with experiment and as shown below that i^*(/> and K* p data cannot be fitted simultaneously by 
two universal parameters pA and (t)A- This may be a potential problem for QCDF. 

The large longitudinal polarization fraction of B K*^ p'^, fi = 0.96to?6' measured by BaBar 
seems to be peculiar as a smaller Jl of order 0.5 is observed in other K*p modes such as K*^p~ 
and K*°p°. At first sight, it appears that the BRY's prediction of fL{K*~p°) = OMtomtoit can 
account for the BaBar measurement. However, as we note in Appendix D, there are sign errors 
in the expressions of the annihilation terms A^'^ and ^3'' (see Eqs. ()D2p and (ID 171) . respectively) 
by BRY: The signs of the r^^ terms in these two equations are erroneous in [13(1 . Because of 
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, I , I , I , I , y t , I , T~~ — 1- — I — I — ^— P 1,1,1, r~~, 1 — r-'i 

0.2 0.4 0.6 0.8 1. 0.2 0.4 0.6 0.8 1. 0.2 0.4 0.6 0.8 1. 

Pa Pa Pa 

(a) (b) (c) 

FIG. 1: Predicted branching ratios [(a) and (b) are for the total and longitudinal branching ra- 
tios, respectively] and (c) the longitudinal polarization fraction for B~ K*^p~ as a function of 
the parameter p^. The solid and dashed curves correspond to the central value and the allowed 
theoretical uncertainty due to the variation of (pA, respectively. The horizontal band represents 
experimental values with one sigma errors. 



the (incorrect) cancelation between r^^ and r^'^ terms in A^'^ , BRY claimed (wrongly) that the 



longitudinal penguin annihilation amplitude is strongly suppressed, while the (3^ term receives 



X X 

^3 is SllUllgiJ' SUppiCSSCU, WllliC UliC fj^ 

sizable penguin annihilation contribution. If a wrong sign for r^^ terms is used, both rates and 
longitudinal polarization fractions will be reduced, especially for the K*^p^ mode where Jl is 
reduced by more than a factor of 2 (see the case (iii) of Table rvT|) . For comparison, BRY's predictions 
are shown in the last two columns of the same table. Using the correct expressions for A^'^ and 
we find that fL{K*-p°) is reduced to the 70% level and /L(i^*V°) is predicted to be 0.39t|]J?. 



The latter agrees with the experimental value fL{K*^p^) = 0.57 it 0.12 64]. As explained above, 
the corresponding prediction 0.221q'q31q'^4 by BRY is too small owing to the incorrect signs in their 
A^^^ amplitudes. 



In short, we have the pattern (see also 13l |) 



fL{K*-p'') > fL{K*-p+) > fL{K*^p-) > fLiK*'p") (4.11) 

for the longitudinal fractions in B ^ K* p decays. Note that the quoted experimental value 
fL{K*~p^) = O.OG^Q ^g in Tables HVl and IVTl was obtained by BaBar in a previous measurement 
where K*-p° and K*- fo{980) were not separated 0]. This has been overcome in a recent BaBar 
measurement, but the resultant value fL{K*^p^) = 0.9ib0.2 has only 2.ba significance j47]. At any 
rate, it would be important to have a refined measurement of longitudinal polarization fraction for 
K*^p^ and K*^p^ and a new measurement of fL{K*~p'^) to test the hierarchy pattern (j4.1ip . 



B K*(t> 

Experimentally, B — > K*<j) decays have been studied with full angular analysis from which 
information on final-state interactions can be extracted. Historically, it was the observation of 
large transverse polarization in these decays that had triggered the theoretical and experimental 
interest in the study of charmless B — > VV decays. 

Theoretically, B K*(j) decays can be analyzed in the same manner as the K*p modes. The 
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decay amplitude of i? K* (f> can be approximated as 



A-.K*^^ ^ + «f + - ^«3,Ew)^]^**- (4.12) 



When the penguin annihilation contribution is turned off, we have 

B-^K*-d> 



I c,0 




(4.13) 



From the amplitude parameters given in Table IVl it is clear that there exists a constructive (destruc- 
tive) interference in the A~ {A^) amplitude. As a consequence, although the factorizable amplitudes 
respect the hierarchy l-'^^*^] : l-'^ft^*,^! : l-'^i^*^! = 1 '■ 0-35 : 0.007 due to the iy — A) structure of 
weak interactions and helicity conservation in strong interactions, the negative- and longitudinal- 
helicity amplitudes are comparable in magnitude. Numerically, we indeed find fi = 0.62 and 
f\\ ~ /± = 0.19 (see Table rvT|) . This is very similar to the decay K*^p^ where fi, is also 

found to be small, of order 0.50 . Experimentally, the naive expectation of /l ~ l — Amy/m'^ ~ 0.90 
is strongly violated in charmless penguin-dominated VV modes. Nevertheless, a small /l for K*(f) 
is quite natural in QCD factorization because the parameters are helicity dependent. The fact 
that real parts of as and 03, ew Aip signs from h = to h = — and that is smaller in magnitude 
for the longitudinal amplitude (see Table |V]) will render the negative helicity amplitude comparable 
to the longitudinal one. 

Even though the longitudinal polarization fraction is reduced to 60% level in the absence of 
penguin annihilation, this does not mean that the polarization anomaly is resolved. As stated 
before, irrespective of the predictions for polarization fractions, the first task we need to focus 
on is to reproduce the correct rate for B — > K*(l) because the calculated branching ratio of order 
4.1 X 10^^ is too small by a factor of ~ 2.5 compared to the measured one, ~ 10 x 10^^ (cf. Table 
lIVp . Assuming weak annihilation to account for the discrepancy between theory and experiment, 
we can fit the data of branching ratios and /l simultaneously by adjusting the parameters pA 
and 4>A- However, this also means that QCDF loses its predictive power in this manner. We find 
that the rate and fi can be accommodated by having pA ~ 0.65 and (pA ~ —53°. This set of 
the annihilation parameters differs slightly from that extracted from B — > K*p decays, namely, 
Pa{K* p) ~ 0.78 and (j)A{K*p) ~ —43°. Therefore, within the framework of QCDF, one cannot 
account for all charmless B VV data by a universal set of pA and (pA parameters. This could 
be an indiction that large penguin annihilation cannot be the ultimate story for understanding 
B VV decays. 

Since the complete angular analysis oi B K*(j) has been performed by both BaBar and 
Belle, information on the parallel and perpendicular polarizations and their phases relative to the 
longitudinal one is available. We see from Table ITVl that f± and /y are very similar, of order 0.25. 
Experimentally, the phases (p\\ and (j)± deviate from either vr or zero by more than 4.6(T and 3.3cr 
for K*^(f) and 5.5cr for K*^(j) This implies the presence of final-state interactions. The relative 
phases are calculated to be 

</.||(K*-0) = (80ii)° (expt: (46 ± 10)°), 

</>y(i?*V) = (781^?)° (expt: (44^?)°), (4.14) 
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and (/)\\{K*- (/)) - cI)_l{K*- c/)) - vr = (j)\\{K*° c/)) - cI)_l{K*°4') - tt k, 0.7°. They are consistent with the 
data. 

Thus far we have chosen the renormahzation scale to be /i = mf,(mb) in calculations. We now 
address the issue with fi. In principle, physics should be independent of the choice of n, but in 
practice there exists some residual n dependence in the truncated calculations. We have checked 
explicitly that the decay rates without annihilation are indeed essentially stable against fi. How- 
ever, when penguin annihilation is turned on, it is sensitive to the choice of the renormahzation 
scale because the penguin annihilation contribution characterized by the parameter 63 is domi- 
nantly proportional to asifJ'h)c6{iJ-h) at the hard-collinear scale fih = V fJ-^h- For the hadronic 
scale A,, « 500 MeV, we have /i/, « 1.45 GeV and 1 GeV for /u = 4.2 GeV and 2.1 GeV, re- 
spectively. At the amplitude level, the enhancement of penguin annihilation at = 2.1 GeV is of 
order Qs(l)c6(l)/[as(1.45)c6(1.45)] ~ 1.8. We find that if the renormahzation scale is chosen to 
he n = mfe(mb)/2 = 2.1 GeV, we cannot fit the branching ratios and polarization fractions simul- 
taneously for both B — > K*(j) and B K* p decays. For example, the rate of the former can be 
accommodated with pA ~ 0.25, but the corresponding //, ~ 0.28 is too small. Likewise, although 
B{B~ — > K*^p~) can be fitted well with pA ~ 0.55, the resultant fi < 0.12 is highly suppressed. 
This is ascribed to the fact that at the scale p = 2.1 GeV, the negative-helicity amplitude receives 
much more enhancement than the longitudinal one and hence the longitudinal polarization is sup- 
pressed at the small p scale. In order to ensure the validity of the penguin-annihilation mechanism 
for describing B VV decays, we will confine ourselves to the renormahzation scale p = mh{miy) 
in the ensuing study. 



B K*uj 



The decay amplitudes for B K*uj read 



V2A 



B- 



X' 



K* 



X 



K* 



V2A 



X 



h _ 



+ 



1 



X 



(4.15) 



Prom the previous analysis of K*p and K*(f) decays, we found two distinct sets of the penguin 
annihilation parameters pA and (pA- If the set of parameters inferred from B K* p decays, 
namely, pA = 0.78 and cjiA = —43°, is employed, we obtain B{B~ —>■ K*~uj) « 4.5 x 10~^ and 
B(K K*^uo) w 3.9 X 10-^ which are slightly higher than the respective experimental upper 
limits, 3.4 x 10~^ and 2.7 x 10^^ [5a, |61[. By contrast, if the parameters pA = 0.65 and (pA = —53° 
extracted from i^*</> modes are used, the resultant predictions B{B^ K*^uj) ~ 3.5 x 10^^ and 
B{B^ — > K*^uj) ~ 3.0 X 10^^ are consistent with experiment (see Table irv|l . Of course, if the B ^ to 
form factors are smaller than what we expected as implied by the measurement of B~ P~oj, 
then B{B K*uj) will be safely below the current bounds. 



B K*K* 
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The expressions of B 
Ab-^k*ok*- 



K*K* decay amplitudes read 



V ^ p=u,c 



Gp 



1 



p=u,c 



v/2 



2 '^4,EW + + /^3,EW 



^{BK\K')^ 



K 



2 "4,EW 



T"4,EW + 03 + 04 ~ o/^3,EW 



p=u,c 



+ fnfh 



^4,EW 



(4.16) 



Both 



^0 



K*^K*^ and K*^K*~ are h ^ d penguin-dominated decays, while -B" 

K*~K*~^ proceeds only through weak annihilation. Hence, their branching ratios are expected 

Recentlj^the K*°K*^ mode was first measured by BaBar with 
Our prediction is slightly smaller, about la away 
from the BaBar measurement (Table ITVl) . The absence of transverse polarization in the K*^K*^ 
mode is due to the approximation we have adapted; that is, we have neglected the transverse 
annihilation contributions relative to other terms. Hence, transverse polarization does not 

^2 ' "4 

Eq. ^M)]- 



to be small, of order < 10 
the branching ratio (1.28lo;32) x 10"^ [59] 



receive contributions from the annihilation terms 6]*^ , 6^ , 6j , ^4^^-^^ within our approximation [see 



Comparison with other works 

Within the framework of QCD factorization, we have studied charmless B — > VV decays closely 
to the works of Kagan ll[ and Beneke, Rohrer and Yang (BRY) IJ]. Nevertheless, there are some 
differences between our work and theirs as we are going to discuss below. 

Without penguin annihilation, Kagan found ~ 0.90 for the K*^(j) mode, while BRY got 
/l ~ 0.67 and we obtained /l ~ 0.62. Kagan did not consider vertex corrections and hard spectator 
interactions in his realistic calculations of a^'^, though he has discussed the latter briefly. It seems 
to us that are essentially helicity independent in the Kagan's calculation and this accounts for 
the difference in the estimation of fi- Moreover, what is the initial value of fi in the absence of 
penguin annihilation is immaterial because we will use the unknown annihilation parameters to 
accommodate the data of branching ratios and fi rather than to predict them. 

We differ from BRY mainly for using different pA_ and (pA_ parameters for describing B — > K*p 
decays. If we follow BRY to use the set of pA and (pA parameters extracted from B — >■ K*(j) decays 
to describe K*p modes, the rates for the latter will be systematically below the measurements. 
Since it is necessary to reproduce the measured rates first in order to have a reliable estimate of 
polarization fractions, we need to fit the K*p data separately. The resultant pA and (j)A parameters 
differ from the ones determined from K*(j) modes. This can be viewed as a potential problem of 
QCDF. 

In the pQCD approach, the calculated branching ratio of B ^ K*(p is too large, of order 15x 10^^ 
and fi ~ 0.75 14|. It has been proposed in 43(] that a smaller form factor Aq^* (0) ~ 0.30 will 



bring down both the rates and fi and bring up /y and f±. While this sounds plausible, the NLO 
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corrections to helicity-dependent coefficients should be taken into account in this approach as 
we have demonstrated that NLO corrections to a~ and will bring down fi significantly. It is also 
important to compute the rates and polarization fractions for B K*p decays in this framework 
and compare with experiment. 

Another plausible solution is to consider the long-distance rescattering contributions from some 
charm intermediate states such as D^*^Di*^ [igI . [itI . 1S|. The idea is simple: First, B D*D* 
decays are CKM favored and hence final-state interactions via charm intermediate states will bring 
up the K*(j) rates. This is welcome since the short-distance predictions of the branching fractions 
for penguin-dominated B VV modes in most of the models under consideration are too small 
compared to experiment. Second, large transverse polarization induced from B D*D* will be 
propagated to cpK* via final-state rescattering. The unknown parameter in the model for final-state 
rescattering is fixed by the measured rate 18|. However, this approach has one drawback. That 
is, while the longitudinal polarization fraction can be reduced significantly from short-distance 
predictions due to final-state interaction effects, no sizable perpendicular polarization f± is found 
owing mainly to the large cancelations occurring in the processes B — > D*D (j)K* and B 
DgD* (f)K* and this can be understood as a consequence of CP and SU(3) symmetry ISj. That 
is, final-state rescattering leads to the suppression of /t(= /± + /||) at the expense of /y » As 
pointed out in 1^, one possibility to circumvent the aforementioned cancelation is to consider the 
contributions from the even-parity charmed meson intermediate states. In view of the fact that 
even at the short-distance level, fi can be as small as 40% to 70% after NLO corrections to effective 
Wilson coefficients are taken into account, it is worth re-examining this type of solution again. 

In soft-collinear effective theory (SCET) [l^, large transverse polarization in penguin-dominated 
VV modes may arise from the long-distance charming penguin contribution. Indeed, the aforemen- 
tioned mechanism of final-state rescattering of charm intermediate states mimics the charming 
penguin in SCET, while both QCDF and pQCD approaches rely on penguin annihilation to resolve 
the polarization anomaly.^ 



B. B ^ VA decays 

The calculated branching ratios and longitudinal polarization frac- 
tions for the decays B Ap, AK*, Alo, A(p with A = 
ai(1260),6i(1235),/s:i(1270),Ki(1400),/i(1285),/i(1420),/ii(1170),/ii(1380) are collected in 
Tables EmXl 



1. B ^ aiV, biV decays 



The decays 5^ (a^,5]*')(p ,7r ) are governed by the decay constants of the p and vr, re- 



spectively. Since fp » f-j^, we thus expect to have B(B^ a^p ) » B{W a'^ir ) and 



Ways of distinguishing penguin annihilation from rescattering have been recently proposed in 63 1 
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B(B^ btp') > B(W bfiT') ~ 10 X 10-^ These features are borne out in our realistic 
calculations (see Table rvTI|) . Calderon, Muiioz and Vera (CMV) [211] found the other way around: 
{af ,bf )p~ modes have rates smaller than {af ,bf)TT~ ones, which we strongly disagree. Since the 
modes (a]~, 6^ )(/)+, vr"*") are governed by fa^ and /fe^, respectively, and since fa^ ~ fp and /^^ is very 
small (vanishing for the neutral we anticipate that a^p"*" and aj'vr"'' have comparable rates and 
the bip^ mode is highly suppressed relative to the b^ p" one. ^ The decays (aj~,6j~)/j° receive both 
color-allowed and color-suppressed contributions: 

oc (a? + . . ^'"^ + (4 + • • ■)xi''^'^'\ 

oc (a? + • • ''^^ + («2 + • • ■)Xr'''\ (4.17) 

where x^^^^'^^^) i^j^e factorizable amplitude defined by Eq. (I3.16p . Since the color-allowed 
amplitude of the b'^ p^ mode is highly suppressed by the smallness of the bi decay constant and 
the color-suppressed amplitude is suppressed by the small ratio of 02/01, it is clear that B{B^ 
CLi p^) ^ B{B~ bip^). The decays (aj~,6|f)tj should have rates similar to {a^ ,bi)p^ . So far 
there is only one experimental measurement of i? — > VA decays, namely, — > with the 



result 



B{B° afp'^)B{at (37r)±) < 61 x 10"^ (4.18) 

Assuming that af decays exclusively to p^ir^, we then have the upper limit of 61 X 10"^ for the 
branching ratio of B^ — > afp^. Our prediction B{B^ af p^) ^ 59 x 10~^ is thus consistent with 
experiment. 

To discuss the effect of the annihilation contribution, let us take the penguin-dominated decays 
bt)K* as an example. From Eq. (A3) of [l| we have 

■^BO^a+K*- °^ («4 + «4,EW)-'^^^"'''^ ^ + ^fBfajK*{bs - ^&3,Ew)aii^* , 

Abo_,+k'- ^ («4 + a2,Ew)^^^'^'''*^ + ifsftjK^bl - hl^^),,K': (4.19) 

where we have replaced the decay constant f^,-^ by f^^ as explained before (see the paragraph after 
Eq. (|2.52p ). From penguin-dominated B — > VV decays we learn that the predicted rates in default 
are typically too small by a factor of 2 ~ 3. In the absence of the experimental information for 
penguin-dominated B VA decays, we shall use the penguin-annihilation parameters pA = 0.65 
and (pA = —53° inferred from B — > i^*(/> decays as a guidance for annihilation enhancement. 
Since the magnitude of 63 is large for the biK* modes (specifically, b'^(biK*) = — 1.78 + 9.92i and 
b^{aiK*) = — 0.19 + 4.11z), B biK* decays receive more enhancement from penguin annihilation 
than B — > aiK* ones. When penguin annihilation is turned off, we have, for example, 

B(^ ^ atK*~) = i3.6t\i1:'oi) x B{^ ^ ft+K*") = {A.ltHtH) x lO'^, 

B{B~ ^ a^i?*0) = {A.ltiltll) X 10-^ B{B- ^ 6^i?*0) = (4.0l|°l°;^) x lO^^. (4.2O) 



As explained in [l|, within the QCD factorization approach, the suppression of 6j~(7r+,p+) modes is not 
directly related to the smallness of the hi decay constant, but is ascribed to the tiny coefficient ai. However, 
the smallness of ai has to do with the decay constant suppression. For more details, see [l|. 
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TABLE VII: Branching ratios (in units of 10~^) and the longitudinal polarization fraction (in 
parentheses) for the decays B — > (ai, &i)(p, (j), K*) with oi = ai(1260) and hi = 6i(1235). 
The theoretical errors correspond to the uncertainties due to variation of (i) Gegenbauer moments, 
decay constants, quark masses, form factors, the parameter for the B meson wave function, and 
(ii) PA,H-, 4>A,H-, respectively. For longitudinal polarization fractions, we add all errors in quadrature 
as the theoretical uncertainty is usually dominated by (ii). Default results are for pA = 0.65 and 
(j)A = —53°. We use the light-cone sum rule results for the B ^ ai and B ^ h\ transition form 
factors (see Table HH) . The model predictions by Calderon, Muhoz and Vera (CMV) [21[ are also 
displayed here for comparison. 
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We see from Table IVlIl that aiK* and biK* modes are substantially enhanced by penguin annihi- 
lation. Experimentally, it is thus very important to measure them to test the importance of the 
penguin annihilation mechanism. 

We have checked explicitly that, in the absence of penguin annihilation, the longitudinal polar- 
ization fractions are close to one half in aiK* modes and 90% in hiK* ones. This can be seen from 
the ratio of the negative- and longitudinal-helicity amplitudes 

(4.21) 



c,0 
1^4 



4,EW 



^Bai,K* 
Bai,K* 



As discussed in the section of i? ^ VV decays, the interference between a^^ and 04 ew generally 
constructive for /i = — and destructive for h = 0. Since \X% „,| : iXg ,,,| : iXt, .^,1 = 1: 
0.50 : 0.06, and \X%. : \X^, : ^1 = 1: 0.21 : 0.03, it is obvious that the A' 

' I Bbi,K*' ' Bbi,K*' ' Bbi,K*' ' 

amplitude of aiK* channels has more chance to be comparable to ^4*^ than the biK* ones. When 
penguin annihilation is turned on, it is evident from Table IVTIl that aiK* modes are dominated by 
transverse polarization amplitudes, whereas biK* are governed by longitudinal polarization states. 
The decays B — > (ai,6i)(/> are highly suppressed relative to the tree-dominated {ai,bi)p modes 
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as they proceed through h ^ d penguin process and are thus suppressed by the small coefficients 
for penguin operators. Moreover, they do not receive any annihilation contribution! 

On the experimental ground, our calculations suggest that the tree-dominated channels af p~ , 
tti p~ , a\p~ , ai p^, afw, bip~ and bip~ should be readily accessible to B factories. Measurements 
of the penguin-dominated modes aiK* and biK* are crucial for testing the mechanism of penguin 
annihilation. 



2. B ^ Ki{1270)V, Ki{im)V decays 

To obtain the branching ratios and fi listed in Table rvTIII for B KiV decays, we have used 
the light-cone sum rule results for the B a and B Kib form factors given in Table IHl The 



decays B — > Ki(f> have been considered in [20|] based on the generalized factorization framework 
where nonfactorizable effects are lumped into N^^, the effective number of colors. It is interesting 
to note that the results of 201] for B Ki{l270)(l) are similar to ours when A'^^^ is close to 5, but 
the predicted rates for Ki{1400)(f) are smaller than ours irrespective of the value of N^^. From Eqs. 
(|B10p and (jBlip we have the decay amplitudes given by 

Ah r c , c I nc-iYiBKi(1270),4,) 

oc [ai + /Jslx,,,^'"'"^^ sin 9k, + [a^ + /^alx,,,^''''^" ^os Oj,, , 
•^B-^K-{im'P i"3 + a4 + P3J^h (4-22) 
OC [ai + al + /JsIk,^,^''''^^ ^os 9k, - [a^ + al + f3'^,^^,F^'''^ sin0^,, 

where F^^^^ denotes generic form factors for the B Kia transition and likewise for F^^^^ . In 
our convention, form factors F^^^^ and F^^^^ have opposite signs (see Table HI]). Since the mixing 
angle 9ki is negative, it follows that the two amplitudes in Eq. (j4.22p contribute constructively 
to B^ Ki{1270)^(j) and destructively to B^ i^i(14OO)~0. Hence, it is naively expected 
that the former has a rate larger than the latter. Indeed, when the penguin annihilation (P^) is 
turned off, we find B{B- i^i(1270)^(^) « 3.2 x lO^^ > BiB^ i^i(1400)-(^) « 3.1 x IQ-"^. 
However, this feature is dramatically changed in the presence of weak annihilation with p^ = 0.65 
and <pA = —53°. Since (33{KiA<i)) and j3^{KiB4>) are opposite in sign, the interference becomes 
destructive in B^ Ki{1270)~ 4> and constructive in B~ — > i^i(1400)~</). This explains why we 
have B{B~ Ki{1270)-4>) < B{B- Ki{U00y(t>) in Table EHll If this relation is not borne 
out by experiment, this will indicate that the parameter pA and hence weak annihilation are small. 

The decays B — > i^i(1270)p have rates larger than that of i? ^ Ki{14:00)p and this can be 
understood as follows. Their decay amplitudes have the expressions, for example, 

« [al + /?3]pA-i^ fK,A sin 9k, + [at + /3f ]^,^^^ fk,^ cos 9k, , 
« [«2 + l3i],K,A fK,A cos 9k, - [al + Pi]^^^^ Jk^^ sin 9k, ■ (4.23) 



^W^K-{liQO)p+ L 4 I ^-3J h 
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TABLE VIII: Same as Table IVIII except for b ^ s penguin-dominated decays (top) B 
Ki{p, K* ,uj) and b ^ d penguin-dominated ones (bottom) B KiK* for two different mixing 
angles = —37° and —58°. Default results are for pA = 0.65 and (pA = —53°. 



1 l£l/^OT7" 




Oki = 


-37° 


Oki = 


-58° 














B 


h 


B 






K^{l270)p+ 


1« 0+9. 8+54.7 
^^"•"-6.8-13.8 


Q r;7+0.39 
"•"^'-0.30 


iQ 1+12.1+47.6 
8.5-14.8 


n 40+0-48 
"•^^-0.36 


■^^ 


:^i(1270)p0 


Q -1+5.3+34.2 
^•-'--3.5- 8.6 


"•^'-'-0.37 


Q 0+6. 0+30.4 

y-o_4 1_ 8 


"•^"-0.30 


B 




-I 7 n+10-7+53.0 
^ ' 7.3-15.3 


n 52+°-^'^ 

"•"-■^-0.36 


r,n 1+12.5+48.5 
^"•-•-- 8.9-15.2 


0.47t°i^ 


B 


■ iff (1270)p0 


Q n+5.0+20.4 
°-^-3.7- 6.1 


"•"3"-0.34 


in q+5. 9+19.1 
1U.0_4 5_ g 7 


0.561°:^^ 




:^?(1270)w 


7 Q+4. 7+24.0 
'•"3-3.1- 7.5 


n r;n-H0.39 
"•"^■^-0.32 


o 9+5.3+21.5 
°-^-3.7- 7.4 


0.48l°i° 


B- 


■ iff (1270)a; 


7 4+4.7+22.6 
'•^-3.2- 6.9 


U.OD„o.22 


o 7+5.4+20.6 
°- ' -3.7- 7.1 


"•"-'^-0.34 




;^?(127O)0 


q fj+l. 7+4.8 
''•O-l. 3-2.9 


n R7-<-0.33 
"•"3' -0.64 


q 9+2.1+5.2 
"J-^-1. 4-2.7 


"•"J-L-O.Sl 


B- 


■ iff (1270)0 


q + 1.9+5.1 

"^•O-l. 5-3.1 


n fi7-l-0-33 
"•"3' -0.64 


q 1+2.2+5.5 
"-••^-l. 5-2.8 


n qi+0-69 
"•"3 ■'--0.37 




iff (1400)p+ 


o f;+2. 8+13.0 
"•"-2.3- 4.3 


"•"-'*-0.23 


cr 7+1.2+17.5 
' -1.0- 4.7 


n O7+0.09 
0-«'-0.43 




;^?(1400)p° 


11 1+3.2+15.5 
-LJ--^-2.8- 5.8 


U.DO_o.2l 


Q Q+2. 1+22.4 
^•°-2.0- 8.0 


n 90+°-°^ 

u.yu_o 28 


B- 


■K°i{im)p- 


in Q+2.6+15.2 
^U-y_3.i_ 5.4 


U.DO_o.26 


7 r+2.06+19.4 
'•^-1.6- 6.0 


o.sstoil 


B- 


■ iff (1400)p0 


q O + l. 3 + 6.4 

"^•O-l. 1-1.9 


n 61+0-32 
"•Oi_0.22 


1 c;+0.7+7.7 
-•-■^-0.3-1.0 


o.7ot°:l^ 




Zi(1400)cj 


c- f;+2. 8+6.8 
"-'•0-2.2-2.5 




1 (.+2.8+9.5 
*-"-l. 7-3.5 


o.9ol°:°8^ 


B 


■ iff (1400)u; 


A cr + 1. 8+6.0 
^•"J-1.4-2.1 




q 1+1.3 + 8.1 
"-••■•--0.8-2.4 


0.87tn^ 




(1400)0 


in /I +7.9+38.3 
J- '-'•'i-5. 1-10.4 


"•^"-0.02 


in 7+7.1+37.69 
-■-"•' -4.6-10.4 


n ^7+0-31 

"•^' -0.22 


B- 


■ iff (1400)0 


11 1+8.5+41.1 
-Li^i.s 4_ii 4 


U-^^-0.09 


11 q+7.5+40.2 
■'■^•"-'-4.9-11.1 


n ^7+0-32 
"•'^' -0.22 




iff (1270)if*+ 


n ni -fO.01+0.03 

'-'•^■'--0.01-0.00 


1.0 


n nn+0-00+0-01 

"•""-0.00-0.00 


1.0 




ifi+(1270)if*- 


n nfi+0-03+1-43 

"•"-""-0.02-0.02 


1.0 


n nft-i-o.02+0.91 

"•""-0.02-0.00 


1.0 




iff (1400)if*+ 


n ns+0-04+0.28 

"•'-'°-0.03-0.00 


1.0 


n nQ-i-O-05+0.30 

"•"^-0.03-0.00 


1.0 




ifi+(1400)if*- 


n nn+0-00+0.20 

"•""-•-0.00-0.00 


1.0 


n nn-i-O-OH-0.69 

"•""-0.00-0.00 


1.0 




;^i(1270)if*° 


r, 4^-^0.25+0.65 
u.^u_Q 26-0.35 


0.86l°:i 


n 9Q+0.19+0.41 
o.^»_0.20-0.07 


0.84t°i^ 




ifiO(1270)K*° 


n nQ+0-04+3.52 

U.U»_Q g3_Q QQ 


0.34l°|l 


n 9t;+0.09+3.53 
o.^'J_0.08-0.00 


0.521°:^^ 


B 


■ if?(1270)if*- 


n nc;-^o.o4+2.2i 

"•"^-0.02-0.00 


n qq+0.58 
"•"J"3-0.06 


n 1 q+0.04+2.00 
"•■•-"3-0.05-0.00 


0.54l°:i 


B- 


■ iff (1270)if*° 


n 1Q+0.13+0.28 
"•-■■^-0.13-0.15 


0.84^°:°^ 


n 1 t;+0.10+0.20 
"•■'^'-'-0.10-0.03 


o.sotn^ 




;^i( 1400) if *° 


r, r|i +0.01+3.41 
"•"■•--0.01-0.00 


0.97^°:°? 


n ns-i-0-05+3.01 

"•"°-0.06-0.01 


n 89+°-°^ 




if?(1400)K*° 


n cri +0.08+1.30 
"•"J ■•--0.11-0.29 


n 71 

"•'-•--0.29 


n qfi+0.08+1.60 
"•"3O-0.01-0.30 


0-85lo;54 


B- 


■ if?(1400)A'*- 


n 9O+0.04+0.53 
"•^"-0.06-0.14 


n 77+0.20 
'-0.34 


n 1 Q+0.04+0.73 
"•■•^^-0.04-0.10 


0.85t°:l^ 


B- 


■ iff (1400)if*° 


n ni -1-0.00+2.00 

"•"■•■-0.01-0.00 


0.93^°:?^ 


n nc,-^o-03+i.77 

"•"^-0.04-0.01 


0.92tH? 



Just as the case for B ifi0 decays, the interference is constructive (destructive) in iff (1270)/?+ 
and destructive (constructive) in iff (1400)/)"'" in the presence (absence) of weak annihilation with 
PA = 0.65 and (pA = —53°. This explains why the rates of the former are larger than the lat- 
ter, especially for 9ki = —58°, in Table IVIIII Hence, measurements of the relative rates of 
ifi(1270)p and ifi(1400)p will enable us to see the role played by the weak annihilation effect. 
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If B{B Ki{1270)p) < B{B ifi(1400)p) is observed, this will hint at the smallness of weak 
annihilation. The reader may notice that the decay modes involving Ki{\27{)) and -fri(1400) al- 
ways have opposite dependence on the mixing angle 6ki- For example, K^[\27Q)p^ gets enhanced 
whereas K^{l'\QQ)p'^ is suppressed when Ok^ is changed from —37° to —58°. 

Decay rates of i? — > i^i(1270)-fC* and Ki{\4Q{))K* are generally small because they proceed 
through b ^ d penguin diagrams and are suppressed by the smallness of the penguin Wilson 
coefficients. Their branching ratios are of order 10^^ — 10^^. The decay modes K*"^ and 
K^K*^ are of particular interest as they are the only AV modes which receive contributions 
solely from weak annihilation. Just as the decay K*^K*^ discussed before, the absence of 

transverse polarization in the K*^ and K^K*^ modes is due to the fact that the annihilation 
terms bf,bf-,bf, 6^^^ vanish under our approximation. 

From Table it is clear that the channels iCf p+, K^p~ for Ki = i^i(1270) and i^i(1400) 
have sizable rates and the experimental search of them would be encouraging. 

3. B ^ fiV, hiV decays 

Results for the decays B (/i, w, i^*, 0) with /i = /i(1285), /i(1420) and hi = 
/ii(1170), /ii (1380) for two distinct sets of the mixing angles Osp^ and 9ip^ are summarized in 
Table HXl Among tree-dominated decays, the channels /ii(1380)p~ for 6ip^ = 0°, /i(1285)/9~ and 
/ii(1170)p^ have branching ratios of order 10~^ as they receive color-allowed tree contributions. 
Many of the penguin-dominated modes e.g. /i(1420)Er* have branching ratios in the range of 
(5 ~ 15) X 10~^. It is of interest to notice that the decays involving /ii(1380) in the final state have 
a sharp dependence of the rates on the mixing angle 9ip^ . 

C. B ^ AA decays 

For the axial vector mesons ai(1260), &i(1235), /i(1285), /i(1420), /ii(1170), /ii(1380) and 
i^i(1270), i^i(1400), there exist many possible B AA two-body decay channels. We will classify 
them into tree- and penguin-dominated decays. The latter involves the strange axial- vector meson 

1. Tree- dominated decays 

The decay amplitudes for some of tree-dominated B — > AA decays are shown in Appendix B. 
Since the decay constant of the bi is either vanishing or very small, it is anticipated that bibi 
channels are highly suppressed relative to aioi. Some of aibi decays are comparable to aiai. We 
find that aia^ and aj~a5 modes have rates larger than the corresponding p^ p^ and p^ p^ ones, 
but a^a? is very similar to p^ p^ . While 05*^0^, a^a^, 0,1^1 a^b^ modes have branching ratios 
of order (20 ~ 40) x 10^^, the other channels are suppressed by the smallness of either /^^ or the 
coefficient 02. 
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TABLE IX: Same as Table IVTIl except for the decays B (/i, i^*, (/>) with fi = 
/i(1285),/i(1420) and hi = /ii (1170), /ii (1380). We use two different sets of mixing angles: (i) 
(93 = 27.9° and Oip^ = 25.2° (in first entry), corresponding to = -37°, and (ii) Oap^ = 53.2°, 
6*1 = 0° (in second entry), corresponding to = —58°. 



Mode 



h 



Mode 



B 



B- ^ /i(1285)p- 



B° /i(1285)w 



B- fi{1285)K* 



B° fi{1285)K*^ 



^ ^ /i(1285)</. 



B- hi{U70)p' 



B^ /ii(1170)p° 



B^ /ii(1170)u; 



B- hi{ll70)K*- 



W hi{1170)K 



B^ /ii(1170)(/> 



*o 



in 0+5. 5+0.5 
^^•^-3. 5-0.4 
o q+5.1+0.4 
°-^-3.2-0.3 
n 9+0.2+0.3 
U-^-0. 1-0.0 

n n+0. 1+0.3 

U-^-0. 1-0.0 
+1.1+2.5 
-0.4-0.2 
+1.0+2.2 
-0.4-0.1 
+8.3+10.8 
-2.8- 2.3 
h3.8+21.4 
-2.2- 4.8 

c: c+7.9+10.1 
0-0_2.7- 1.7 

cr 1 +3.6+20.0 
'^•■'--2.1- 4.7 

+0.002+0.010 
0.001-0.000 
+0.002+0.009 
0.001-0.000 
+10.1+2.8 
6.6-3.3 
+6.5+1.8 
4.2-2.1 
+0.11+1.16 
0.03-0.00 
Q Q^+0.09+0.78 



0.91 



+0, 



1.0 
0.9 
5.8 
5.7^ 



0.002 
0.002 
17.4 
10.9 
0.05 



0.90l[! 
0.77l[! 
0.711° 

+0 
-0 
+0 
-0 
+0 
-0 
+0 



0.87 
0.86 
0.90 
0.47 



0.89l[! 
0.451° 

+0 



-0.02-0.00 

-I r+2.4+1.2 

-■-•^-0.9-0.2 
+1.4+0.8 
-0.6-0.1 
+2.5+12.8 
-1.6- 4.3 
+5.1+31.6 
-3.0- 7.1 
+2.2+11.5 
-1.4- 4.2 
7 1+5.1+30.1 
'•-■^-2.9- 6.9 



0.9 
5.3 
7.7 
4.5 



0.006: 

0.001 



[-0.007+0.010 
-0.002-0.005 
1-0.005+0.074 
-0.003-0.000 



0.93 
0.90 
0.96 
0.96 
0.22 



-0.68 



+0 
-0 
+0 
-0 
+0 

0.811° 

hO 
-0 



0.99 
0.99 
0.84 



0.82] 
0.811° 



0.97 
0.93 



B- 



b" 



B' 



b" 



B' 



b" 



B- 



./i(1420)p- 
/i(1420)pO 
/i(1420)w 

> fi{U20)K*- 



/i(1420)K*° 



/i(1420)(/. 

> /ii(1380)/9- 
/ii(1380)p° 
/ii(1380)a; 

> /ii(1380)K* 
hi{1380)K 
/ii(1380)(^ 



*o 



n 1+0.2+0.1 

'-'■-'^-0.1-0.0 

1 q+O.6+0.2 

-^•"-•-0.3-0.0 

n ni +0.03+0.02 

-L-0. 00-0. 00 

n n/1+0.12+0.08 

'-'•'-'^-0.03-0.00 



0.02 
0.1 



+0.03+0.05 
-0.01-0.00 
+0.2+0.3 



15.9 
15.6 



-0.1-0.0 
+8.4+18.0 
-5.3- 7.0 
+10.9+10.4 
5.2- 4.7 



0.70 



0.87 
0.53 
0.86 
0.50 
0.64 



1/1 c+8. 0+17.4 

-lyl Q+10. 2+10.1 
5.0- 4.6 



0.0001 

0.0008 



+0.0001+0.0006 
-0.0000-0.0000 
-0.0009+0.0009 



0.9 



-0.0001-0.0001 
+0.5+0.2 



5.9 
0.01 



-0.3-0.2 
+3.2+1.0 
-2.1-1.1 
+0.01+0.04 



-0.00-0.00 
Q 02+0.03+0.40 



0.1 
0.5 
8.1 
3.7 
8.3 
3.9 
0.004 
0.007 



0.01-0.00 
+0.1+0.1 
0.1-0.0 
+0.7+0.4 
0.3-0.1 
+4.0+21.3 
2.8- 6.6 
+2.0+7.8 
1.3-2.2 
4.4+21.8 
-2.9- 6.9 
+1.9+8.3 
1.3-2.6 
+0.003+0.230 
0.002-0.000 
0.005+0.170 
0.003-0.001 



r0.19 

-0.07 

0.931°:°^ 
0.38l°:i 

0.08 
0.40 
0.31 
-0.39 
+0.04 
0.76 
+0.48 
0.52 
+0.37 
0.61 

n 40+0.49 

U-O'i-0.61 
+0.03 
0.75 
+0.02 
0.44 
+0.00 
0.08 
+0.01 
0.07 
+0.59 
0.13 
0.99 
0.00 



0.97 
0.98 
0.95 
0.95 
0.30 
0.01 
0.97 
0.98 
0.87 



+0.02 

-0.14 
+0.01 
-0.12 
+0.13 
-0.75 



0.88l°i^ 
0.88l°i^ 



0.88 
1.00 
0.99 



0.12 

-0.64 
0.00 
0.04 
0.01 
0.14 



Among various B {ai,bi){fi, hi) decays, we see from Table 1X1 that only a^/i(1285) and 
a^/ii(1170) modes and a^/ii(1380) with 6ip_^ = 0° can have sizable rates and all other charged and 
neutral channels are suppressed. 



2. Penguin-dominated decays 

The penguin-dominated B — > AA decays involve at least one Ki meson. Results for the decay 
modes Ki{ai,bi, fi,hi,Ki) are summarized in Table IXll Some salient features are (i) T{B 
Ki(1270)ai) > r{B E:i(1400)6i) > r{B Ki{l270)bi) > r{B E:i(1400)ai), (ii) B 
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TABLE X: Branching ratios (in units of 10 ^) and the longitudinal polarization fraction for tree- 
dominated B — AA decays. For decays involving fi and hi states, we use two different sets of 
mixing angles: (i) 6*3 = 27.9° and Oip^ = 25.2° (in first entry) and (ii) 63 = 53.2°, 61 = 0° (in 

sccoud ou(ry). 



Mode 



h 



Mode 



B 



fL 



B- 
B- 
B- 
B- 



07 0+23.9+11.4 

^'■°-15.3- 5.3 

1 n+1-6+6.2 
-•-•^-0.5-0.1 
1 /1+2. 5+2.8 
-■■•^-1.0-0.0 

q + 5.6 + 11.0 

•J-^-2.3- 0.8 
-1+1.6+15.7 
0.7- 0.3 



0.73t°:i 



1.0^ 



0.95 
0.95 
0.96 



0.00 

-0.82 
0.02 
0.80 
0.03 
0.65 









07 4+16.1+9.7 
-4-13.7-1.4 






a?a? 


n t-+0.8+9.3 
"••J-0.2-0.0 


b" 




a^bf 


41 q+20.7+16.6 
^■•■■^-18.2- 3.4 






0-1 K 


n Q+l.09+3.6 

U.O_g 4_Q ]^ 


b" 




a?6? 


q Q+6.2 + 2.6 

^•"-2.3-0.5 



b" 



0.60ir.g 

^•^'^-0.05 
hO.OO 
-0.80 



0.98 



0.981°:°} 



B- 
B- 
B- 
B' 
B- 
B- 
B- 
B- 



► ar/i(1285) 
. ar/i(1420) 
^ a^/ii(1170) 
^ a^/ii(1380) 
.6r/i(1285) 
. 6r/i(1420) 

► 6^/11(1170) 

► 6i-/ii(1380) 
/i(1285)/i(1285) 
/i(1285)/i(1420) 
/i(1420)/i(1420) 
/i(1285)/ii(1170) 
/i(1420)/ii(1170) 



+5.6+6.9 
4.3- 0.7 



4.1-0.8 
+0.2+0.3 



12.4 

11 o+^-''+^-° 

0-2_o. 1-0.0 
1 c+o. 4+0.9 
-^•^-0.3-0.0 

99 ^+14.5+5.3 
^^•4- 9.3-3.1 
14 1+9.5+3.4 
-'-4- 1-6.0-1.9 
1 9+0.7+0.3 
-^•^-0.5-0.1 



n 7q+0.22 
"3-0.32 

71+0-23 
•J- '1-0.31 



0.42; 

0.77" 



-0.42 
-0.19 
-0.16 
-0.33 



7.6 
0.8 



4.7+1.6 
3.0-0.7 
+1.3+4.3 
0.5-0.3 



n 7+1.0+3.5 
' -0.4-0.1 

n nq+O.17+0.16 
"-'■^<5-0.02-0.01 

n 9+0.6+0.8 
u-^-o.i-o.i 
1 9+2.0+9.2 

-■-•^-0.9- 0.5 
r, 0+1.3+3.0 
^•°-0.3-0.1 

n nc^+o-io+o-34 

'-'•^^-0.04-0.00 

0.5+3.0 
-0.2-0.2 
0.3+3.1 
0.1-0.0 
+0.2+2.5 
-0.1-0.0 
+0.01+0.10 
-0.01-0.00 
0.05+0.63 
-0.02-0.00 
0.00+0.03 
0.01-0.00 



0.82 
0.79 
0.73 



0.3 

0.3 

0.2 
0.01 
0.05 
0.01 

r, r,i +0.01+0.06 
U-UJ^-0.00-0.00 
1 1+1.9+1.1 
-'-•-'--0.7-0.2 

0.6 
0.02 



+1.1+0.6 
0.4-0.1 

+0.06+0.03 
-0.01-0.00 



rvqi+0.02 

"-'•^-'--0.22 

qo+° °2 
u.yu_Q 24 

89+°-°^ 
u.oy_Q 24 

+0.16 
-0.56 
+0.19 
-0.57 
+0.23 
-0.33 

0.89t°:°? 

h0.03 
-0.76 

0.94l°:?i 

o.85t°:°^ 

0.03 
-0.77 
0.06 
-0.84 
+0.07 
0.84 
+0.31 
0.28 
0.10 
0.66 
+0.05 
0.34 

0.68l°:i 
0.981°:°? 

0.02 
0.09 
+0.11 
0.23 



0.95: 



0.94 
0.67 
0.66 
0.26 
0.57 
0.94 



0.97 
0.87 



a?/i(1285) 

a?/i(1420) 

a?/ii(1170) 

o?/n(1380) 

b?/i(1285) 

50/1(1420) 

6?/ti(1170) 

6?/ii(1380) 

/ii(1170)/ii(1170) 

/ii(1170)/ii(1380) 

/ii ( 1380) /ii (1380) 

/i (1285)^1 (1380) 

/i (1420)^1 (1380) 



n 1+0.1+3.1 
"-'•-'--0.1-0.0 



0.1 



+0.1+2.7 



■--0.0-0.0 

n no+0.02+0.12 
^•^^-0.01-0.01 
n no+O-02+0.12 
"-••^^-o.oi-o.oi 
rv 1 +0.3+2.1 

^^•-'--0.1-0.0 

n ns+O-17+1.36 

^•^°-0.06-0.02 

n n-i +0.01+0.12 
"-'•"-'-'--0.01-0.00 
n nf;+o.o7+o.79 

"-'•"-''-'-0.03-0.00 

n 9+0.4+2.7 
"-'•^-0.1-0.1 
n 9+0.3+2.5 
"-'•^-0.1-0.0 
n ni +0.08+0.10 

"-'•^-'--0.01-0.01 

n 1+0.3+0.5 
"-'•-'^-0.1-0.0 
n 9+0.2+5.1 
"-'•^-0.1-0.0 



0.1 



n ni +0.01+0.20 
^•^-'--0.01-0.00 



0.04 



+0.04+1.73 



-0.02-0.00 
Q e+1.6+1.6 



0.4 
0.1 
0.3 



n ns+o-i'''+o.io 

^•U°-0.04-0.01 



0.961°:°^ 



^-0.6-0.3 
+0.7+1.2 
0.3-0.1 
+0.1+0.1 
-0.0-0.0 
0.6+0.7 
0.3-0.1 
n ni +0.01+0.42 

^•^-'--0.01-0.00 

n n«+o.i4+o.79 

^•^°-0.05-0.02 
n r,c+0.09+0.07 
^•U^-0.03-0.01 

+0.6+0.3 
-0.2-0.0 

n nni +0.001+0.003 
u.uui_Q 000-0.000 

o f^4+0. 10+0.08 
^•^^-0.03-0.01 



0.3 



^s+°-i^ 

^•^•3-0.59 



0.49: 



-0.06 
-0.52 

0.80 



0.181 



0.15 



24+°-^^ 
^•^4_o 26 



0.30 



+0.70 
-0.34 



+0.2+3.4 
-0.1-0.0 



0.32l°:^I 

0^08l°:i 
+0.53 

0.38 

U.OD_o 22 
0.29 
-0.49 

0.811°:!^ 

0.861°:^^ 

+0.12 
-0.79 



0.48 



0.66 



0.86 



n 47+0.45 

+0.17 
-0.70 
0.01 
-0.55 
+0.02 
0.60 
+0.01 
0.74 
+0.01 
0.70 
+0.03 
0.39 

96+° °^ 
"-'•^'-'-0.62 

n Q7+0 02 
^•^'-0.22 

0.97l°:°? 
74+0-22 

'"1-0.25 
+0.04 
0.12 



0.81 
0.97 
0.97 
0.97 
0.96 
0.97 



0.95 
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TABLE XI: Branching ratios (in units of 10 ®) and the longitudinal polarization fraction for 
penguin-dominated B — KiA decays for the mixing angles: (i) 9ki = —37°, Osp^ = 27.9° and 
6*1 = 25.2°, and (ii) Ori = -58°, Osp^ = 53.2° and Oip^ = 0°. Default results are for pA = 0.65 

and OA = —53°. 



Decay 



= -37° 



= -58° 









B 


II 


B 


II 






Ki{1270)ai 


AC) 0+58.2+165.6 
^z.o_27 3_ 41.1 


24+°-^^ 

U.Z4_Q Q7 


Aa 1+60.9+176.8 
^u.i_29 8_ 43.7 


l6+°-^^ 

U.lO_o.o6 






(1270)0? 


91 /.+29.7+81.2 
zi.u_^3 g_2o.9 


n 97+0.69 
•J-^' -0.21 


99 /I +30.8+88.7 
'^•^-15.0-22.2 


^•1^-0.03 


B 




■ Kf{1270)a^ 


A A O+60. 9+165.9 
^^•'^-28.6- 43.5 


n 90+0.09 

"•'^"-'-0.19 


10 q+63.4+175.7 
^°-'-*-31.1- 47.0 


0.151°:?^ 


B- 




■ i^f (1270)a? 


90 O+30.1+84.7 
^o-0_14.3_2l.4 


o.26l°:i 


9fi q+31. 6+87.6 
^"•"^-15.6-23.2 


0.201°:^} 








19 n+10.6+12.0 
6.1- 7.9 


0.36l°:i 


7 C-+7.3+21.7 
'•^-3.3- 6.1 


0.96tni 




-> 


i??(1400)a5 


f- 7+5.5+7.4 
"•'-3.3-4.5 


'-'•^'3_o.4i 


c- t;+4.2+12.4 
^•^-2.3- 4.1 


0.98l°ii 


B- 




■ (1400)aj; 


iq + 11.4+17.1 

^o.o_ 6.8-10.1 


42+°-58 
"•^^-0.40 


Q n+8-4+23.8 
»-0_3 9_ 7 4 


0.98lgii 


B- 




(1400)a? 


f. n+5.4+7.0 
"^•"-3.0-4.3 


n QQ+0.35 
U.OO_o.3o 


Q 1+3.6+10.2 
''•J^-1.5- 2.6 


0.941^:°^ 






(1270)6^ 


14 0+12.9+65.8 
i4.0_ 7 2-13.7 


0.281^:^2 




0.131^:11 






^?(1270)6? 


7 0+6.3+33.0 
'•'-*-3.5- 6.7 


""^^-0.18 


r- q+7.4+34.1 
"•^-3.1- 5.2 


"•-•-^-0.16 


B- 




-^?(1270)6r 


1 c- 0+13.8+72.6 
i-'J 'J- 7.5-14.9 


o.3ilg:?J 


10 n+15.1+69.5 
-'-•J-'J- 6.0-10.8 


"J-"JO-0.09 


B 




i^f (1270)6? 


Q n+8.3+36.3 
^•U-4.5- 7.7 


0.39l°i^ 


1+9.2+32.2 
o^-L-S.O- .x9 


n 99+0.65 
'J-^'^-0.22 






i^f (1400)6^ 


9c: n+27.3+205.5 
Z,O.U_ii i_ 22.9 


0.9lini 


9/^ 9+24.2+209.2 
^O.Z_i2.0- 24.7 


^•^^-0.56 






(1400)6? 


10 n+13.4+110.1 
^^■'J- 5.6- 12.3 


n qi+0.05 
"•^-"--0.66 


1 7+12.3+109.4 
-"-■J-' - 6.2- 13.1 


0.981^1 


B 




ifi°(1400)6f 


97 7+28.6+231.4 
^'•'-12.0- 26.1 


n 01 +0-05 

U-y-'-_0.66 


on /I+27.2+235.5 
"^"•^-13.9- 29.3 


0.981°:°? 


B- 




iiTf (1400)6? 


1 O+14.4+107.9 
5.9- 12.2 


92+°-^^ 
"•^^-0.35 


14 4+13.4+112.7 
^-^•^- 6.7- 13.8 


0.991°:?! 






X0(1270)/i(1285) 


14 c;+20.1+68.4 
8.2-11.1 


U.OD_Q 7Q 


c- 7+8.8+31.0 
"^•'-3.7- 5.2 


0.221°:^^ 






K?(1270)/i(1420) 


10.4^1:11^11 




IQ Q+20.9+82.1 
^^"•^- 8.0-10.7 


0.641°:^? 






K?(1270)/ii(1170) 


r q+8.7+25.7 
^•^-3.4- 4.4 


"•^^-0.60 


4 1+11.3+7.8 
^•-1-- 3.4-2.3 


0.83lo:83 






K?(1270)/ii(1380) 


(-+13.3+38.4 
5.6- 5.2 


0.931°:?^ 


cr+11.5+33.5 
"•^- 4.9- 4.3 


0.501°:^ 


B- 




iff (1270)/i (1285) 


1 c; 7+21.5+73.4 
-L<J- ' _ 8.7-11.7 


0.601°:^^ 


P_ 9+9.0+34.2 
u.^„3 8- 5.2 


29+°-^° 
u-z»_0 35 


B- 




i^f(1270)/i(1420) 


in Q+9-7+53.7 

iu.y_3 7_ g 4 


n qQ+0.09 
'-'•^"^-0.55 


IQ 7+21.8+85.6 
^^•'- 8.2-11.0 


0.651°:^^ 


B- 




(1270)/ii(1170) 


R r-+9.3+27.4 


U.OD_Q 52 


cr ^+12.2+9.2 
4.0-0.8 


0.85l°:^i 


B 




ivTf (1270)/ii(1380) 


Q 9+14.7+40.4 
6.1- 5.8 


u.yo_Q 73 


Q+12. 3+33.5 
"•^- 5.2- 4.4 


U.OU_o 73 






K0(1400)/i(1285) 


A n+5. 3+12.2 
^•U-2.1- 4.5 


'-'•-'-^-0.13 


9 9+2.2+10.7 
^•^-0.5- 1.7 


"•^^-0.56 






X? ( 1400) /i (1420) 


10 Q+25.4+87.6 
J-O-y_16.7-l9.0 


0.041°:^! 


91 c+26.2+95.4 


0.261°:!^ 






K?(1400)/ii(1170) 


Q Q+10.4+75.4 
4.5- 8.3 


0.961°:?^ 


91 n+21-6+148.3 
^-■-•"-11.4- 19.7 


o.9olg:°? 






( 1400) /il (1380) 


Ifi K+22-0+81.9 
-■■"•"-10.1-16.9 


0.55l°:l? 


(? c;+8.4+33.2 
"•"J-3.9- 6.5 


"•^^-0.08 


B- 




E:r(1400)/i(1285) 


A q + 5.9+13.2 

^•o_2.4_ 4.4 


"•-■-^-0.10 


9 9+2.5+10.4 
^•^-0.5- 1.7 


U.OU_o 73 


B- 




i^f ( 1400) /i (1420) 


IQ /I +26.2+89.7 
-■-^•^-13.0-19.8 


0.03lHt 


99 C-+27.3+99.9 
^^•0-14.1-18.1 


^•^'-0.15 


B-^ 




i^f (1400)/?,i(1170) 


q 1+11.1+74.5 
4.6- 8.3 


0.96ini 


99 n+23.2+154.5 
^^•"-12.1- 20.8 


n ni +0.08 
"•^-■--0.83 


B- 




ivTf (1400)/ii(1380) 


1 7 7+23.4+87.8 
•'-10.8-18.0 




7 n+9.1+35.9 
/.U_4 2_ 6.9 


^•^^-0.18 
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Table XL (Continued) 



Decay 




Ok, = -37° 




0K, = -58° 


B 


fL 


B 


fL 



B- 
B- 
B- 
B 



(1270)Kf (1270) 
-ftTf (1270)^1+ (1400) 
i^f (1400)^1+ (1270) 

(1400)^1+ (1400) 
X?(1270)i^'?(1270) 
K?(1400)A'?(1400) 
K?(1270)if?(1400) 
ii:?(1270)i??(1400) 

> E:?(1270)Kf (1400) 
^ -fCf (1270)i^5'(1270) 

> -ftTf (1270)^^10(1400) 

> i^f (1400)i^?(1400) 



0.07tg; 



0.01 
0.20 



0.041° 
0.161°: 
0.29lg 
0.85lg; 

o.oolg 

0.03l|] 
13+°- 
0.79l[j 



0.21 



+0 
-0. 



12+1.07 
05-0.00 
02+0.35 

01- 0.00 
20+3.72 

11- 0.00 
05+0.42 

02- 0.00 
17+1.03 

12- 0.01 
23+4.99 
21-0.01 
44+10.95 
41- 0.38 
01+19.0 

00- 0.0 
02+17.1 

01- 0.0 
11+0.43 
07-0.02 
42+8.01 
39-0.29 
16+3.16 
12-0.02 



1 
1 

0.671°:^ 
0.64l°:i 
0.77l°:l 

"■""-0.70 

-0.03 
-0.86 
+0.23 
-0.40 

0.74l°;l 



0.94: 
0.51 



0.10 
0.01 
0.23 



0, 
0. 
+0, 
-0, 
+0, 
0, 

0.04l°: 
0.44l°: 

o.oelg 

0.541°: 

o.o2i°: 

-0, 
-0, 
+0, 



0.41 



+0.28 
-0.47 



o.io: 

0.23 
0.46 
0.11 



10+1.23 

06- 0.00 
02+0.34 
01-0.00 
24+3.03 
17-0.00 
08+0.63 
03-0.00 
33+8.35 
31-0.00 
06+0.13 
01-0.0 
37+4.17 
27-0.23 
02+14.12 
01- 0.0 
05+13.07 

07- 0.00 
17+4.28 
13-0.02 
30+2.29 
21-0.13 
10+0.36 
05-0.04 



1 
1 

"•'"-0.73 

0.461°:^^ 
0.83l°:^i 

"•^"-0.38 

+0.04 
-0.91 
+0.26 
-0.51 

r, 70+0.20 

U- '9-0.84 
+0.34 
-0.49 



0.91 

0.42 



0.61 



Ki{ai,bi) decays are dominated by transverse polarization amplitudes except for i^i(1400)6i and 
Ki(1400)ai with Ok^ = —58°, and (iii) the charged and neutral B decays have similar rates and 
longitudinal polarization fractions. For example, B{B~ — > {fi,hi)) « B{B^ Ki{fi,hi)), 
B{B- i^f (a?,6;)) B(B° K^{a\,hl)) and B{B- i^O(af ,6^)) ^ B{^ i^f (a+ 6+)). 

The first feature can be understood as follows. Consider the decays B° — (al,6l) as an 
illustration. Their decay amplitudes are given by 

1 



Ah 

'■^-^ft:-(1400)o+ 
Ah 

^K^ {I270)b+ 

4<, 



oc 

oc 
cx 

(X 

oc 



+ a 



c,h 
4,EW 



c,h 



)/^3,EW 



aiKi 



h 



[•••], J Ki^/i^iA sinOK, + [•••], 
[■ ■ -L^K.^fKiA COS Ok, -[•••] 

[■ ■ ■i.K.^fK.A sin^i^i + [• 

[■■■Lk,Jk,a cos9k^ - [■ 



,K,sfK,B COS 9k„ 



■>i^j7(i400)fe+ ' JbiKi^^"iA "1 L JOiKjs-'^is "i- (4.24) 

Since Kitti modes are dominated by transverse amplitudes and since the negative- helicity param- 
eters such as a'^ {uiKia) and (aiKiB) have opposite signs, it is clear that the interference is 
constructive in B ^ ^1(1270)01 and destructive in B — *■ iCi(1400)ai for a negative mixing angle 
6ki ■ This also explains why the former increases and the latter decreases when the mixing angle 
is changed from —37° to —58°. For the Ki(1400)6i modes dominated by the longitudinal ampli- 
tudes, they have large rates as a^^ibiKiA) and a^{biKiB) are of the same sign. In general, Kiai 
and Kibi rates are insensitive to the value of Okh —37° or —58°, except for ifi(1400)ai modes. 
It is interesting to notice that i<'i(1400)ai channels are dominated by transverse amplitudes for 
9ki = —37° and by longitudinal ones for = —58°. Therefore, measurement of polarization 
fractions in i? ^ iiri(1400)ai will yield a clear discrimination between the two different Kia — Kib 
mixing angles. Branching ratios of -B — s- Ki{fi,hi) fall into the range of 10~^ ~ 10~^. At first 
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sight, it appears that they depend on the mixing angles and 9ip-^ or Bsp-^. However, the latter 
two angles are correlated to the first one [see Eq. (|2.4|) ]. Consequently, the decays B — > Ki{fi, hi) 
depend on only one mixing angle Oki- From Table [Xll it is clear that the mode -ftTi ( 1400) /ii (1380) 
has a strong dependence on 6ki ■ 

Just as B ^ KiK* decays, branching ratios of the h ^ d penguin-dominated KiKi modes are 
small, of order 10^'' — 10^^ owing to the smallness of the penguin coefficients. 

In short, there are many penguin-dominated B A A decays within the reach of 

B factories: A'i(1270)ai, A'i(1400)6i, A:i(1270)6f, Ki{im)af , Ari(1270)(/i(1285), /i(1420)) 
and i^i(1400)(/i(1420), /ii(1170)). In most cases, transverse polarization is large except for 
A"i(1400)(6i,/ii(1170)), A"i(1270)(/i(1420),/ii(1380)) with Ok^ = -37° and A"i(1400)ai with 
= —58° where longitudinal polarization dominates. 

V. CONCLUSIONS 

In this work we have presented a detailed study of charmless two-body B decays into final 
states involving two vector mesons (VV) or two axial-vector mesons {AA) or one vector and one 
axial- vector meson iVA), within the framework of QCD factorization, where A is either a ^Pi or 
^Pi axial-vector meson. Owing to the G-parity, the chiral-even two-parton light-cone distribution 
amplitudes of the ^Pi (^Pi) mesons are symmetric (antisymmetric) under the exchange of quark 
and anti-quark momentum fractions in the SU(3) limit. For chiral-odd light-cone distribution 
amplitudes, it is other way around. The main results are as follows. 

• We have worked out the hard spectator scattering and annihilation contributions to P ^ VA 
and B A A decays. 

• NLO nonfactorizable corrections to longitudinal- and negative-helicity effective Wilson coef- 
ficients generally differ in magnitude and even in sign. For some VV modes, the construc- 
tive (destructive) interference in the negative-helicity (longitudinal-helicity) amplitude of the 
B VV decay will render the former comparable to the latter and bring up the transverse 
polarization. Any serious solution to the polarization puzzle should take into account NLO 
effects on a^. 

• The measured rates and /l of penguin-dominated charmless VV modes can be accommodated 
(but cannot be predicted at first place) in QCD factorization by allowing for sizable penguin 
annihilation contributions. However, the parameters pA and (f>A fit to the data of K*(/) and 
K*p are not the same. Hence, we do not have a good hint at the values of pA and (pA for 
B AV and B ^ AA decays. 

• While NLO contributions due to vertex, penguin and hard scattering corrections are insen- 
sitive to the choice of the renormalization scale p, the penguin annihilation contribution at 
the hard-collinear scale is sensitive to p. In the present work, we choose p = mi,(mi,) for 
the reason that if = mh{mh)/2 is selected, the decay rates and polarization fractions of 
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B K*4> OT B ^ K*p cannot be simultaneously fitted by the annihilation parameters pA 
and (j)A- 

The predicted rates and longitudinal polarization fractions by QCD factorization for tree- 
dominated pp modes are in good agreement with experiment, but the calculated B{B~ — ^ 
p~io) is slightly high. The latter may imply that the B ^ u transition form factors are 
slightly smaller than what are expected from the light-cone sum rules. Only in the decay 
B^ p^LO where a large deviation from the naive expectation of ~ 1 is possible. We 
found ~ 0.55. 

Using the measured K*^p~ channel as an input, we predict the branching ratios and polar- 
ization fractions for other B K* p decays and find the relation fL{K*~ p^) > fL{K*^p^) > 
fL{K*^p^) > fL{K*^p^). Experimentally, it is quite important to measure them to test 
theory. Our result of fL{K*^p^) 0.39 is consistent with experiment and is higher than the 
prediction ~ 0.22 made by Beneke, Rohrer and Yang. 

The calculations suggest that the tree-dominated channels a,f p~ , a^p~, a^p~ , cliP^, aiUJ, 
bfp~ and b^p^ should be readily accessible to B factories. One of the salient features of 
the ^Pi axial vector meson is that its axial- vector decay constant is small, vanishing in the 
SU(3) limit. This can be tested by measuring various bip modes to see if TiB" ^ b^p+) < 
r(S° bfp-) and T{B- b]; p^) < T{B- b\p-). 

In the absence of the experimental guideline, we employed the penguin annihilation param- 
eters PA = 0.65 and (pA = —53° inferred from the channel B K*(j) to describe penguin- 
dominated B VA, AA decays. It is very crucial to measure the penguin-dominated modes 
aiK* and biK* to test the importance of penguin annihilation. We found that aiK* modes 
are dominated by transverse polarization amplitudes, whereas biK* are governed by longi- 
tudinal polarization states. 

For B KiV decays involving the Ki{\27Q) or iiri(1400) meson, the channels p^ , K^p~ 
have sizable rates and the experimental search of them would be encouraging. Measurements 
of the relative strengths of Ki{127G)(f){p) and Ki(lA{)G)(f){p) will enable us to test the impor- 
tance of weak annihilation. The rates oi B ^ ^1(1270)^* and Ki{lAOQ)K* are generally 
very small. The decay modes K*^ and K^K*^ arc of particular interest as they are the 
only VA modes which receive contributions solely from weak annihilation. 

Among the decays B {fi,hi){p,u, K* with /i = /i(1285), /i(1420) and hi = 
/ii(1170),^i(1380), the tree-dominated modes hi{1380)p', fi{1285)p~, hi{1170)p- and sev- 
eral of the penguin-dominated channels e.g. /i(1420)K* have appreciable rates. 

For tree-dominated B — AA decays, the afa^, a^Ui, a^b^ and a^b^ modes have sizable 
branching ratios, of order (20 ~ 40) x 10~^. Among various B — (ai, 6i)(/i, hi) decays, only 
a]~/i(1285) and aj'/ii(1170) modes and a]^/ii(1380) with Oip^ = 0° can have large rates and 
all other charged and neutral channels are suppressed. 
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• There are two salient features for penguin-dominated B A A decays: (i) T(B 
Ki(1270)ai) > r{B Ki{im)bi) > T{B i^i(1270)6i) > r{B Ki(1400)ai) and 
(ii) most of them are dominated by transverse polarization amplitudes except for i^i(1400)6i 
and -fri(1400)ai with 9ki = —58°. Since the Xi(1400)ai channels arc dominated by trans- 
verse amplitudes for O^i = —37° and by longitudinal ones for = —58°, measurement 
of polarization fractions in B ^ Ki{l400)ai will yield a clear discrimination between the 
two different Kia— Kib mixing angles. Many penguin-dominated B AA decays are 
are readily detectable at B factories: ^1(1270)01, ifi(1400)6i, Ki(1270)6f , ifi(1400)a^, 
Ki(1270)(/i(1285),/i(1420)) and Ki(1400)(/i(1420), /ii(1170)). 
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APPENDIX A: FLAVOR OPERATORS 



The coefficients of the flavor operators a^'^ can be expressed m terms of a^'^ (sj, [fil as fohows:^ 

a^(MiM2) = a^(MiM2), 
a^{MiM2) = 4(MiM2), 

h,Pf A^^^MiM^) - 4'P{MiM2) for M1M2 = VA, AA, 



al'^{MiM2) + al'^{MiM2) for M1M2 = VV, AV, 



^h,p,^ M 1 = i < (^1^2) + 2 a^'^(MiM2) for M1M2 = AF, VA, 

^ ^ \ a^'P(MiM2) - rf 2 a^'P(MiM2) for M1M2 = , FF, ^ ' 



a^|w(MiM2) 



f a9'P(MiM2) - a7'^(MiM2) for M1M2 = FA, AA, 
\ ag'P{MiM2) + a)'^{MiM2) for M1M2 = VV, AV, 

h,P M \ = \ ''ioiMiM2) + rf 2 a^'f(MiM2) for M1M2 = AF, VA, 
"^'^w^ ' \ a^5^'(MiM2) - rf 2 a^.P(MiM2) for M1M2 = AA ,VV. 

Note that the order of the arguments of {M1M2) and (M1M2) is relevant. For vector mesons 
we have 

rU,) = ^IM^ (A2) 



while for axial-vector mesons we have 



-» = ^4^. (A3) 

^ 'mbW Ja 

APPENDIX B: DECAY AMPLITUDES 

For simplicity, here we do not explicitly show the arguments, M\ and M2, of oP/^ and 
coefficients. The order of the arguments of of (M1M2) and /3f (M1M2) is consistent with the order 
of the arguments of X^^*^!'*^^)^ where 

/3f (M1M2) = 'J^lMllMlifi, (Bl) 

The decay amplitudes for (ai, 6i)p, (ai, (/i, (/i, hi)K* , Kip and KiK* can be obtained 

from Appendix A of [1] by the replacement of P by y with the same quark content. 

For VA modes, we only list those channels involving uj and^ vector mesons. For other B VA 
decay amplitudes, the reader is referred to Appendix A of [1] with a simple replacement of the 
pseudoscalar meson by the vector meson. 



B ^a^uj /2 P 

V ^ p=u,c 



Spu {a^2 + P2) + 2al''^ + af'^ + ^a|;^w 



The numerical values of the coefficients ai{MiM2) also depend on the nature of the initial-state B meson. 
This dependence is not indicated explicitly in our notation. The same remark applies to the annihilation 
coefficients 5f defined below. 
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9 "4,EW ^ ^^3 ^ 



EW 



+ 



Gf 



X 



{Bu},ai) 



4,EW 



'^^3 2^-'*'EW 2^4,EW 



-^(Bai,a;) 



+ 



Opu{ «2 Pi;t^"4 r,"3,EW 9'^4,E 



:,EW 



2^3,EW 2^"^'^^ 



X. 



{Bu),a\) 



A 



Gf 
V2 



p=u,c 



1 



"3 ~ 2^3'EW 



for B — ai(a;, 
2^0 



p,h 



2"3,EW 



(Bai,0)| 
^(Bai,</))| 



= ^ 



p=u,c 



^"4.EW "I" A'S "T 



1 



+ 



2"4,i!.w ■ /-d ' -^'4 2^3,'ew + 2'^'l'EW 
<^p.(a^ + Pi) + 2al'' + c4'' + ^a|;^w - J«!:ew + P^ + ^ 



1 p,h 



oP,h 



,p,h 



^ RP''' -L ^ RP'h 1 yiBojJ!) ^f7;\r>,h _ ^ p,h 1 y{Buj,ff) 
9/^3,EW "I" 9'^4,EW ^fe "I" V ^ "3 9"3,EW "^/i 



a 



p,h 



p=u,c 



2"3,EW 



+ VBfnU\pl^-hl%^ 



for B ^ 

V2^1 



(B2) 



(B3) 
(B4) 
(B5) 



(B6) 



(B7) 



~Jk 1^ \ [OpuKoti + j + 04 + a4^EW + P3 + P3,EWJ ^/i 
V ^ p=u,c L 

1 

2' 



+ [<5p„at + 2af''^ + k;y4^''^''^4' 



(B8) 



V2^ -0 



/9 ) ["4 2 4,EW"i"''^3 2'^3,EWj"^/i 

V p=u,c K 



+ 



[6pua^ + 2c^/ + ^(4%^ 



X, 



{BKi,u>) 



(B9) 



Ah 

•^b--^k: 



^ 2^ ^P i ["puP2+(^3 +"4 ~ 9O!: 



1 ,P,/» 

2 --3,EW 
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J- p,h 



X 



(BIO) 



^/2 ^ 



p,h , p,h 

tto +04 — —a 



'4 2"3,EW 2""^'^^ 



"•"^3 2'^3,EW 



(Bll) 



for B Ki(u;,(j)). 

The relevant decay amplitudes for B — > AA decays are 

p,h I 3 

2"3,EW 2"4,EW 



» p=u,c 



-t^(Sai,ai) 



(B12) 



V p=u,c 



"r^P4 2'^3,EW "1" 2^4,EW 



-^^(Boi,ai) 



(B13) 



-•^|o_o,o = % E aS.-^) L (at - - a^-'^ + ^c^;^w + - /3f - 2/3f 



2'"4,EW 



-^(Bai,ai) 



(B14) 



V ^ p=u,C \ L 



2 "4,EW P3,EW 



p,ti ^ p,h 

dr, — -OCa 



+ 



,EW 



2 "4,EW 



X 



(B15) 



l/3P,/l 1 ^(-^-^1 '"l) _L ./O 

+P3 ~ 2^3,EW ~ 2^^'EWj^/i 



"3 2 4,EW 



■^.(SaiJf) 



+ 



3 

2'"3,Ji;W 2 



2 "3,EW 2 "4,EW 
•9 



,/oP,/i 1 x-C-BaiJ?) I 

"rP3 2'^3,EW 2'^'^'EWj^/i f' 



(B16) 



for B ai(ai, /i), and 



% E ( \^pM + ) + + <:ew + /3f + /^I'ew 

V ^ p=u,c I 



^(Bai,Ki) 



+ 



p,h 



2 "3,EW 



^(Bifi.ai) 
ft 



11 V-^ p=u,c 



,EW 



^(Bai,i^i) 



(B17) 
(B18) 
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'p=u,c 
V ^ p=u,c 



qP,h 



p,h 



p,h 



f^'" I Z^)'" oPjh I ^ /Oj 



y{Bai,Ki) 



(B19) 



Gf 



X 



{BKi,ai) 
h 



(B20) 



V p=u,c ( 



p,h aP,h np,h 1 ^{BKtJl) I 



(B21) 



V p=u,c 



a 



p,h 



1 

-a 



4,EW "I" o/^3,EW 



2 4,IiVV ' ' d 2 

2 



X 



+ [V"2 + + 2"3,EwJ + V 2 



p.h I p,/i 

03 + a 



p,h 



4 2"3,EW 



2 "4,EW '^3 



1 oP-.h 
2'^3,EW 



h 



(B22) 



for i? — > ii'i(ai,/i), where Ap'^^ = VpbV*^, X^^ = VpbV*s, and the hehcity-dependent factorizable 

,(BMi,M2) 



ampUtudes x^'"^''"^^ are defined in Eq. (j3.16p . The decay ampUtudes for B — > bi{u>,(j)) are 
obtained from B — > ai(a;,(^) by replacing ai ^ 61. Likewise, the expressions for B — > hi{uj,(p) 
decay amphtudes are obtained by setting (/icj hiuj) and (/i0 ^ /iii?i>). 



APPENDIX C: AN EXAMPLE OF THE ANNIHILATION AMPLITUDES IN 
B ^ AV DECAYS 






(c) 



(d) 



FIG. 2: Annihilation contributions to the decay amplitude of B ^ AV , where (a) and (b) corre- 
spond to (c) and (d) give rise to A\^. 



In this appendix we show an explicit evaluation of the annihilation diagrams in Fig. [2] with 
(Mi,M2) = {A^V) and the conventions py ~ En^ and ~ En^. The longitudinal annihilation 
amplitudes of Figs. W^a) and[2)[|b) read 
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^Fig.2(a) = -2{0\d{l-j^)b\B){igsf^^ 



X 1^ dxdyi-l)Tr[M^f{y)^sM^\ Wla{l + 75) 



i{-i){k + pvY 



{k +pv) {xpv + ky 



(CI) 



k=ypA 



and 



A 



Fig 



2(b) = -2(0|d(l -75)6|^)fe)'^^ 



1 rl 



X / / a!xdy(-l)TV| -7„75M|f(y)7'^M|f(x)(l + 75) 



(/c + ypA^ik + paY 



■ (C2) 



The longitudinal projectors and Mil' are given in Eqs. (|2.28|) and (|2.3U|) . respectively. 
Case 1: Taking 

Jv mv{e*v ■ n+) , , 



M\{x) 



4 2 



^iV (x) 



2E 



2 o-^i'Ts (y) - 75 — ^ — 



and using 



(0| J(l - 75)6|^) = i 



fBm% 
mb+md 



(C3) 



(C4) 



we have 



A 



(1) 

Figs.2(a)+2(b) 



Cp f f 

-iT^as-^ fBfvfAr^ / / dxdy^Y{x) 
jVc ^ JO JO " 

-2i7Tas^fBfvfAri [' (\xdy^\{x)^a{y) f4 

-'Vc Jo JO \X 



+ (y) 



(C5) 



where use of Jq dy ^a{y) = has been made. 
Case 2: Taking 



MYix) 



Jymy myje* ■ n+) j i ^ ^ u) 



2E 



h!l{\x) 



and 



we obtain 



A 



(2) 

Figs.2(a)+2(b) 



^11 (y) ^ -i-^ 



(C6) 



(C7) 



iT^as— fsfyfAr; 



h^)^x) + ^hp{x) hY>{x) 



10 JO 
"1 rl 



dxdy $11 (y) 



+ 



Case 3: If 



2i7ras^fBfvfAr^J^ dxdy <i>f{y)'^,{x) (^^^ 
Jy my{e* ■ n 



M\\ {x) 



fi-^Ux), 



(C8) 
(C9) 
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and 



then we have 



A^^^ - -2iTTa — f 

^Figs.2(a)+2{b) - J- 



c 



Case 4: For 
and 

we have 

<.2(b) = 2^^as^JBfvfArl fj\xdy^f{y)^^{x)[^^^ 
Finally, we obtain 

and 



We are led to 



4\AV) = vra, [\xdy \ -r^^a{x)^\{y)'^^^fi- 
J 



+r^$f(x)$.(y)^il±^ 

A II ^Zy 



Likewise, 



2(1 + y) 



xy2 



A II 



-r^^<(x)$„2(y)^^^|, 
A II 2;^y I 

Ai'\VA) = nas £dxdylr];Mx)<^f{y)^^^^^ 



xy^ 



This implies 



(jVV _ (jVA _ _(jAV _ _(jAA _ ^ 
jjVV ^ jjAA ^ _j^AV ^ _j^VA ^ 

which lead to the third line of Eq. (|3.39p . 



50 



APPENDIX D: EXPLICIT EXPRESSIONS OF ANNIHILATION AMPLITUDES 



The general expressions of the hehcity-dependent annihilation amplitudes are given in Eqs. 
(j3.29p - (j3.38p . They can be further simplified by considering the asymptotic distribution amplitudes 
for ^v,^v, ^l^^,^a^^ and the leading contributions to ^l^^,^^^^: 



^{u) = 6uu, ^^^^'{u) = 6uu, $[|^^(n) = 18ai'''''nn(2n - 1), 
^>„(ii) = 3(2^ - 1), <^l^'{u) =3a^'^^{6u'^ -6u + l), <^]f^{u)=3i2u-l), 
<^>]_{u) = 6uu, '^T'i^) = I8a^'^^uu{2u - 1), '^Y'i^) = ^^u, 
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1 ^f{v) 
dv — 
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(D33) 
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ISvra, 



ISvraq 



'Pi 



+ 4 - vr^ 



(D44) 



(D45) 



^i,0^1 D_ 3 ; 



4,0/3 
2,1 1 



A 'Pi] 



(D46) 



where the logarithmic divergences occurred in weak annihilation are described by the variable Xa 

r^^^X,, f'^^-W (D47) 

JO u Jo u 2 

Following these variables are parameterized as 

XA = ln(^]{l + PAe"^--), (D48) 



with the unknown real parameters pA and (pA- For simplicity, we shall assume in practical calcu- 
lations that are helicity independent, = X^ = X^. 



Note that while our result for ^^"(FF) is in agreement with Kagan and Beneke et al. [l 
(up to a sign), the relative sign between r^^ and r^^ in aI'^{VV) {A^^^ (VV)) is positive (negative) 
in our case [see Eqs. ()D2p and ()D17p ] and in but negative (positive) in jl3l |. 



APPENDIX E: EXPLICIT EXPRESSIONS FOR HARD SPECTATOR TERMS 

Using the asymptotic distribution amplitudes, the explicit expressions of the integrals Jq dudv 
appearing in the transverse hard spectator interaction amplitudes H~ {M1M2) and Hf {M1M2) [see 
Eqs. (|3.9p - (j3.13p ] are summarized in Table IXlIl 
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TABLE XII: The explicit expressions for the integrals /q dudv • • • appearing in the transverse hard 
spectator interaction amplitudes H~ {M1M2) and Hf {M1M2) described by Eqs. (j3.9|) - (j3.13p . where 



a = ' \ /? = ^ and the upper (lower) sign is for Hf {H^). 

M1M2 

V1V2 9 

V^Pi 9(X--1) 

V 27a(XJj - 2) 9 

V^Pi -3/3{Xjj-l) 9 -3/3 

^PiV 9{Xjj-l) 3/3 ^3/3 

^Pi^Pi -9a/3(X^-2) 9 -3/3 

'Pi'Pi 9(X^-1) ^3/3 

27a{Xj^ - 2) 

iPi^Pi -3/3(X^-l) 3/3 -6/32 
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